i=1
for(i in 1:length(dependentlist)){
  dep=dependentlist[i]
  print(dependentlist_eng[i])
  load(paste("../data/modeldata/RForigmodeldatawithgeoandgeom_",dep,".RData",sep=""))  
  preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/RF_geoandgeom/RFwithgeoandgeom_fw_5fold_6p_",dep,"_allpreds",sep=""),kk=1:5,endround = 6,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
  print(preds)
  importance_ranfor_pset(modeldata=origmodeldata,dependent=dep,pset=5,altdata=origmodeldata)
  predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
  predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=5,kappasum = F,tausum = F,altdata=origmodeldata)
  predictors = vector()
  for (n in 1:ncol(preds)){
    for (n2 in 1:nrow(preds)){
      predictors <- c(predictors,as.character(preds[n2,n]))
    }
  } 
    print("###################### WITH  PREDICTORS from the FW SELECTION ###################")
  uniquepredictors <- unique(predictors)
  predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=5,kappasum = F,tausum = F,altdata=origmodeldata)
  
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"
## [1] "Prediction error at end is:  0.683982683982684"
## [2] "Prediction error at end is:  0.60995670995671" 
## [3] "Prediction error at end is:  0.637662337662338"
## [4] "Prediction error at end is:  0.628571428571429"
## [5] "Prediction error at end is:  0.61948051948052" 
## [6] "Prediction error at end is:  0.60995670995671" 
##                                                    k 1
## 1                                     geom_10m_fl4_L80
## 2                                      profc_ws7_hr_hr
## 3                                      geom_10m_fl8_L8
## 4        geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr
## 5                              vectorruggedness_hr_ws5
## 6 geom_hr_L50m_fl10_rplippatchdensity_UE_hr_40cells_hr
##                                             k 2
## 1                                           SGU
## 2 geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr
## 3                                 slope_ws7_10m
## 4                              geom_10m_fl3_L22
## 5                              geom_10m_fl1_L20
## 6                              geom_10m_fl2_L20
##                                            k 3
## 1                             geom_10m_fl3_L80
## 2                              geom_10m_fl3_L3
## 3            sagaTopographic_Wetness_Index_50m
## 4                           dtm_hr_TEXTURE_r10
## 5                              crosc_ws3_hr_hr
## 6 geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr
##                                                 k 4
## 1                                dtm_hr_TEXTURE_r10
## 2                                  geom_10m_fl3_L50
## 3 geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr
## 4                                  geom_10m_fl1_L19
## 5                                       TRI_hr_ws52
## 6                                  geom_10m_fl1_L25
##                                               k 5
## 1                                 geom_10m_fl3_L7
## 2 geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr
## 3    geom_hr_L50m_fl1_rplipshape_UE_hr_40cells_hr
## 4                           Profile_Curvature_50m
## 5      geom_hr_L3_fl1_rplipshape_UE_hr_60cells_hr
## 6                                  maxic_ws11_10m
## [1] "OBB error with all predictors of  allpreds is  0.611111111111111"
##                                   MeanDecreaseGini MeanDecreaseAccuracy
## TRI_hr_ws31                              0.2823200          0.001872578
## vectorruggedness_hr_ws15                 0.3234946          0.001773886
## sagaTopographic_Wetness_Index_50m        0.4709617          0.001348061
## TRI_hr_ws30                              0.2515826          0.001338887
## geom_10m_fl4_L9                          0.3161175          0.001287507
## Slope_10m                                0.1602523          0.001238871
## TRI_hr_ws18                              0.2284034          0.001194503
## TRI_hr_ws14                              0.2118763          0.001185284
## vectorruggedness_hr_ws59                 0.2145666          0.001178964
## TRI_hr_ws35                              0.1562579          0.001171159
##         
## altpreds  1  2  3  4  5
##        1  4  0  0  0  0
##        2  0 13  0  0  0
##        3  0  0 21  0  0
##        4  0  0  0 38  0
##        5  0  0  0  0 32
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.583333333333333"
## [1] "OOB-error:  0.527777777777778  for predictors geom_10m_fl4_L80 AND profc_ws7_hr_hr"
##      
## preds  1  2  3  4  5
##     1  2  0  0  0  0
##     2  1  9  0  1  2
##     3  0  0 10  0  0
##     4  1  3 11 36  8
##     5  0  1  0  1 22
## [1] "Kappa overall =  0.62114430869723"
## [1] "Tau overall =  0.664351851851852"
## [1] "mean quality =  0.5485649997704"
## [1] "The quality is  0.5485649997704"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.6794198163761"
## [1] "classification error rate with altdata:  0.268518518518519"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.592592592592593"
## [1] "OOB-error:  0.546296296296296  for predictors geom_10m_fl4_L80 AND profc_ws7_hr_hr AND geom_10m_fl8_L8 AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND vectorruggedness_hr_ws5 AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_40cells_hr AND SGU AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND slope_ws7_10m AND geom_10m_fl3_L22 AND geom_10m_fl1_L20 AND geom_10m_fl2_L20 AND geom_10m_fl3_L80 AND geom_10m_fl3_L3 AND sagaTopographic_Wetness_Index_50m AND dtm_hr_TEXTURE_r10 AND crosc_ws3_hr_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr AND dtm_hr_TEXTURE_r10 AND geom_10m_fl3_L50 AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr AND geom_10m_fl1_L19 AND TRI_hr_ws52 AND geom_10m_fl1_L25 AND geom_10m_fl3_L7 AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipshape_UE_hr_40cells_hr AND Profile_Curvature_50m AND geom_hr_L3_fl1_rplipshape_UE_hr_60cells_hr AND maxic_ws11_10m"
##      
## preds  1  2  3  4  5
##     1  4  0  0  0  0
##     2  0 13  0  0  0
##     3  0  0 21  0  0
##     4  0  0  0 38  0
##     5  0  0  0  0 32
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"
## [1] "Prediction error at end is:  0.711255411255411"
## [2] "Prediction error at end is:  0.722077922077922"
## [3] "Prediction error at end is:  0.731168831168831"
## [4] "Prediction error at end is:  0.686147186147186"
## [5] "Prediction error at end is:  0.676190476190476"
## [6] "Prediction error at end is:  0.695238095238095"
##                          k 1
## 1    geom_dtm_10m_hyd_fl5_L6
## 2   vectorruggedness_hr_ws19
## 3    vectorruggedness_hr_ws7
## 4 terraintexture_hr_ws33_tp5
## 5   geom_dtm_10m_hyd_fl5_L37
## 6            geom_10m_fl8_L6
##                                                  k 2
## 1                          geom_dtm_10m_hyd_fl5_L100
## 2     geom_hr_L3_fl10_rpliprichness_UE_hr_20cells_hr
## 3       SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 4                                    minic_ws5_hr_hr
## 5 geom_hr_L50m_fl10_rplipedgedensity_UE_hr_5cells_hr
## 6                                  geom_10m_fl10_L23
##                                              k 3
## 1   geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr
## 2                    terraintexture_hr_ws33_tp25
## 3   geom_hr_L3_fl1_rplipsimpson_UE_hr_60cells_hr
## 4   geom_hr_L50m_fl10_rplipshape_UE_hr_5cells_hr
## 5     geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr
## 6 geom_hr_L50m_fl1_rpliprichness_UE_hr_5cells_hr
##                           k 4
## 1            geom_10m_fl8_L15
## 2  CrossSectionalCurvature_hr
## 3 Relative_Slope_Position_50m
## 4            geom_10m_fl2_L90
## 5  profc_DTM_50m_avg_ws11_50m
## 6           geom_10m_fl2_L100
##                                                  k 5
## 1                                 dtm_hr_TEXTURE_r10
## 2                        terraintexture_hr_ws45_tp25
## 3 geom_hr_L50m_fl1_rplippatchdensity_UE_hr_5cells_hr
## 4                                   geom_10m_fl4_L14
## 5   geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr
## 6      geom_hr_L50m_fl1_rplipshannon_UE_hr_5cells_hr
## [1] "OBB error with all predictors of  allpreds is  0.583333333333333"
##                                                MeanDecreaseGini
## vectorruggedness_hr_ws49                              0.1989343
## vectorruggedness_hr_ws27                              0.2927074
## Slope_10m                                             0.1633677
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1        0.3949265
## geom_10m_fl4_L9.1                                     0.2475790
## vectorruggedness_hr_ws25                              0.2640938
## TRI_hr_ws21                                           0.2033363
## terraintexture_hr_ws29_tp25                           0.1932838
## maxic_ws7_hr_hr                                       0.2178423
## vectorruggedness_hr_ws9                               0.2809349
##                                                MeanDecreaseAccuracy
## vectorruggedness_hr_ws49                                0.001724241
## vectorruggedness_hr_ws27                                0.001463012
## Slope_10m                                               0.001343616
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1          0.001274637
## geom_10m_fl4_L9.1                                       0.001228478
## vectorruggedness_hr_ws25                                0.001228461
## TRI_hr_ws21                                             0.001201919
## terraintexture_hr_ws29_tp25                             0.001170891
## maxic_ws7_hr_hr                                         0.001121009
## vectorruggedness_hr_ws9                                 0.001067817
##         
## altpreds  2  3  4  5
##        2 33  0  0  0
##        3  0 23  0  0
##        4  0  0 35  0
##        5  0  0  0 17
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.518518518518518"
## [1] "OOB-error:  0.472222222222222  for predictors geom_dtm_10m_hyd_fl5_L6 AND vectorruggedness_hr_ws19"
##      
## preds  2  3  4  5
##     2 28  6  9  1
##     3  2 15  1  0
##     4  3  1 25  3
##     5  0  1  0 13
## [1] "Kappa overall =  0.654502369668246"
## [1] "Tau overall =  0.666666666666667"
## [1] "mean quality =  0.616452991452991"
## [1] "The quality is  0.616452991452991"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.702116508572032"
## [1] "classification error rate with altdata:  0.25"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.537037037037037"
## [1] "OOB-error:  0.462962962962963  for predictors geom_dtm_10m_hyd_fl5_L6 AND vectorruggedness_hr_ws19 AND vectorruggedness_hr_ws7 AND terraintexture_hr_ws33_tp5 AND geom_dtm_10m_hyd_fl5_L37 AND geom_10m_fl8_L6 AND geom_dtm_10m_hyd_fl5_L100 AND geom_hr_L3_fl10_rpliprichness_UE_hr_20cells_hr AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND minic_ws5_hr_hr AND geom_hr_L50m_fl10_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl10_L23 AND geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr AND terraintexture_hr_ws33_tp25 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_60cells_hr AND geom_hr_L50m_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_5cells_hr AND geom_10m_fl8_L15 AND CrossSectionalCurvature_hr AND Relative_Slope_Position_50m AND geom_10m_fl2_L90 AND profc_DTM_50m_avg_ws11_50m AND geom_10m_fl2_L100 AND dtm_hr_TEXTURE_r10 AND terraintexture_hr_ws45_tp25 AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_5cells_hr AND geom_10m_fl4_L14 AND geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_5cells_hr"
##      
## preds  2  3  4  5
##     2 33  0  0  0
##     3  0 23  0  0
##     4  0  0 35  0
##     5  0  0  0 17
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1a.3 - Habitat for soil organisms"

## [1] "Prediction error at end is:  0.462337662337662"
## [2] "Prediction error at end is:  0.407359307359307"
## [3] "Prediction error at end is:  0.397835497835498"
## [4] "Prediction error at end is:  0.435497835497835"
## [5] "Prediction error at end is:  0.426406926406926"
## [6] "Prediction error at end is:  0.37012987012987" 
##                                                k 1
## 1                                      TRI_hr_ws14
## 2                             NormalizedHeight_10m
## 3                                 slope_ws11_hr_hr
## 4                                       TRI_hr_ws5
## 5 geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr
## 6                        vectorstrength_hr_ws11_hr
##                                              k 2
## 1                                            SGU
## 2                      profc_DTM_50m_avg_ws3_50m
## 3                                geom_10m_fl2_L5
## 4                             dtm_hr_TEXTURE_r10
## 5   SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr
##                             k 3
## 1                           SGU
## 2              geom_10m_fl3_L33
## 3     slope_DTM_50m_avg_ws5_50m
## 4 Channel_Network_Base_Level_hr
## 5       geom_dtm_10m_hyd_fl5_L5
## 6               geom_10m_fl3_L6
##                                              k 4
## 1                                            SGU
## 2                      minic_DTM_50m_avg_ws9_50m
## 3                               geom_10m_fl4_L40
## 4                               geom_10m_fl3_L33
## 5                       vectorruggedness_hr_ws23
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr
##                                                    k 5
## 1                                                  SGU
## 2                                        elev_ws15_10m
## 3 geom_hr_L50m_fl10_rplippatchdensity_UE_hr_10cells_hr
## 4          geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr
## 5         geom_hr_L3_fl1_rpliprichness_UE_hr_5cells_hr
## 6                           terraintexture_hr_ws29_tp5
## [1] "OBB error with all predictors of  allpreds is  0.37962962962963"
##                             MeanDecreaseGini MeanDecreaseAccuracy
## elev_ws15_10m                      0.8549385          0.006913781
## elev_ws3_10m                       0.9015790          0.005765188
## elev_ws5_10m                       0.6681540          0.004935603
## Normalized_Height_50m              0.6191737          0.004833930
## vectorruggedness_hr_ws7            0.6394537          0.004193813
## elev_ws11_10m                      0.5455804          0.003863798
## ChannelNetworkBaseLevel_10m        0.4663922          0.003795904
## SGU                                0.6421372          0.003768861
## vectorruggedness_hr_ws27           0.5064353          0.003511857
## NormalizedHeight_10m               0.5839938          0.003378291
##         
## altpreds  1  2  3  4
##        1 24  0  0  0
##        2  0 34  0  0
##        3  0  0 26  0
##        4  0  0  0 24
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.388888888888889"
## [1] "OOB-error:  0.435185185185185  for predictors TRI_hr_ws14 AND NormalizedHeight_10m"
##      
## preds  1  2  3  4
##     1 24  0  0  0
##     2  0 34  0  0
##     3  0  0 26  0
##     4  0  0  0 24
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.416666666666667"
## [1] "OOB-error:  0.259259259259259  for predictors TRI_hr_ws14 AND NormalizedHeight_10m AND slope_ws11_hr_hr AND TRI_hr_ws5 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr AND vectorstrength_hr_ws11_hr AND SGU AND profc_DTM_50m_avg_ws3_50m AND geom_10m_fl2_L5 AND dtm_hr_TEXTURE_r10 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr AND SGU AND geom_10m_fl3_L33 AND slope_DTM_50m_avg_ws5_50m AND Channel_Network_Base_Level_hr AND geom_dtm_10m_hyd_fl5_L5 AND geom_10m_fl3_L6 AND SGU AND minic_DTM_50m_avg_ws9_50m AND geom_10m_fl4_L40 AND geom_10m_fl3_L33 AND vectorruggedness_hr_ws23 AND geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr AND SGU AND elev_ws15_10m AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_10cells_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_5cells_hr AND terraintexture_hr_ws29_tp5"
##      
## preds  1  2  3  4
##     1 24  0  0  0
##     2  0 34  0  0
##     3  0  0 26  0
##     4  0  0  0 24
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1a.4 - Habitat for crops"
## [1] "Prediction error at end is:  0.177056277056277"
## [2] "Prediction error at end is:  0.195670995670996"
## [3] "Prediction error at end is:  0.222943722943723"
## [4] "Prediction error at end is:  0.240692640692641"
## [5] "Prediction error at end is:  0.232034632034632"
## [6] "Prediction error at end is:  0.232034632034632"
##                                                k 1
## 1                                 geom_10m_fl10_L5
## 2                                  slope_ws3_hr_hr
## 3                                geom_10m_fl4_L140
## 4       geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr
## 5   geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## 6 geom_hr_L3_fl10_rplipedgedensity_UE_hr_5cells_hr
##                                          k 2                         k 3
## 1                           geom_10m_fl10_L6            geom_10m_fl10_L6
## 2                           fischerk_ws27_hr             slope_ws5_hr_hr
## 3                           fischerk_ws29_hr   crosc_DTM_50m_avg_ws5_50m
## 4                  planc_DTM_50m_avg_ws3_50m               elev_ws11_10m
## 5                          PlanCurvature_10m             crosc_ws3_hr_hr
## 6 geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr CrossSectionalCurvature_50m
##                  k 4                           k 5
## 1   geom_10m_fl10_L6              slope_ws29_hr_hr
## 2    slope_ws3_hr_hr               slope_ws5_hr_hr
## 3   geom_10m_fl8_L23 DiurnalAnisotropicHeating_10m
## 4 geom_10m_fl10_L110         Mid_Slope_Positon_50m
## 5 geom_10m_fl10_L140              minic_ws29_hr_hr
## 6   geom_10m_fl3_L32                    TRI_hr_ws8
## [1] "OBB error with all predictors of  allpreds is  0.157407407407407"
##                    MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws3_hr_hr           1.2920548          0.011096913
## Slope_hr                  1.2351728          0.010711403
## slope_ws5_hr_hr           1.1428811          0.008930896
## TRI_hr_ws1                1.0350235          0.008450871
## MRVBF_hr                  0.8705280          0.007082092
## TRI_hr_ws2                0.7726378          0.006599331
## TRI_hr_ws5                0.8846360          0.005547737
## geom_10m_fl10_L6.1        0.6122312          0.004701244
## slope_ws5_10m             0.6775568          0.004596555
## slope_ws19_hr_hr          0.5462017          0.004592360
##         
## altpreds  3  4  5
##        3  8  0  0
##        4  0 62  0
##        5  0  0 38
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.175925925925926"
## [1] "OOB-error:  0.111111111111111  for predictors geom_10m_fl10_L5 AND slope_ws3_hr_hr"
##      
## preds  3  4  5
##     3  5  0  0
##     4  2 61  5
##     5  1  1 33
## [1] "Kappa overall =  0.840078973346495"
## [1] "Tau overall =  0.875"
## [1] "mean quality =  0.778019323671498"
## [1] "The quality is  0.778019323671498"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.820806690747015"
## [1] "classification error rate with altdata:  0.0833333333333333"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.175925925925926"
## [1] "OOB-error:  0.12962962962963  for predictors geom_10m_fl10_L5 AND slope_ws3_hr_hr AND geom_10m_fl4_L140 AND geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl10_L6 AND fischerk_ws27_hr AND fischerk_ws29_hr AND planc_DTM_50m_avg_ws3_50m AND PlanCurvature_10m AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND geom_10m_fl10_L6 AND slope_ws5_hr_hr AND crosc_DTM_50m_avg_ws5_50m AND elev_ws11_10m AND crosc_ws3_hr_hr AND CrossSectionalCurvature_50m AND geom_10m_fl10_L6 AND slope_ws3_hr_hr AND geom_10m_fl8_L23 AND geom_10m_fl10_L110 AND geom_10m_fl10_L140 AND geom_10m_fl3_L32 AND slope_ws29_hr_hr AND slope_ws5_hr_hr AND DiurnalAnisotropicHeating_10m AND Mid_Slope_Positon_50m AND minic_ws29_hr_hr AND TRI_hr_ws8"
##      
## preds  3  4  5
##     3  8  0  0
##     4  0 62  0
##     5  0  0 38
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.1 - Average precipitation retention capacity"
## [1] "Prediction error at end is:  0.667099567099567"
## [2] "Prediction error at end is:  0.665367965367965"
## [3] "Prediction error at end is:  0.693506493506494"
## [4] "Prediction error at end is:  0.759307359307359"
## [5] "Prediction error at end is:  0.741125541125541"
## [6] "Prediction error at end is:  0.694372294372294"
##                                    k 1
## 1                      geom_10m_fl4_L3
## 2                          TRI_hr_ws35
## 3                     geom_10m_fl10_L8
## 4                    geom_10m_fl10_L14
## 5                    geom_10m_fl10_L12
## 6 VerticalDistancetoChannelNetwork_10m
##                                        k 2                           k 3
## 1                 vectorruggedness_hr_ws55       geom_dtm_10m_hyd_fl5_L6
## 2 Vertical_Distance_to_Channel_Network_50m Channel_Network_Base_Level_hr
## 3                 vectorruggedness_hr_ws27              dtm_hr_CONVEX_r5
## 4                            elev_ws15_10m       geom_dtm_10m_hyd_fl5_L3
## 5                         geom_10m_fl10_L6                 Convexity_10m
## 6            DiurnalAnisotropicHeating_10m     crosc_DTM_50m_avg_ws7_50m
##                          k 4
## 1           crosc_ws23_hr_hr
## 2            geom_10m_fl2_L4
## 3 terraintexture_hr_ws53_tp5
## 4            slope_ws5_hr_hr
## 5          Plan_Curvature_hr
## 6             minic_ws11_10m
##                                              k 5
## 1                     terraintexture_hr_ws37_tp5
## 2                     terraintexture_hr_ws41_tp5
## 3                  DiurnalAnisotropicHeating_10m
## 4                    terraintexture_hr_ws29_tp25
## 5                     Relative_Slope_Position_hr
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## [1] "OBB error with all predictors of  allpreds is  0.62962962962963"
##                                                MeanDecreaseGini
## MaximalCurvature_10m                                 0.27979971
## vectorruggedness_hr_ws49                             0.24700344
## vectorruggedness_hr_ws59                             0.33433980
## geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr       0.35498139
## geom_10m_fl3_L100                                    0.08331568
## vectorruggedness_hr_ws37                             0.19041487
## geom_10m_fl4_L140.1                                  0.17079245
## vectorruggedness_hr_ws57                             0.18353024
## TRI_hr_ws47                                          0.24291640
## vectorruggedness_hr_ws35                             0.18113058
##                                                MeanDecreaseAccuracy
## MaximalCurvature_10m                                   0.0011014796
## vectorruggedness_hr_ws49                               0.0010956129
## vectorruggedness_hr_ws59                               0.0010255557
## geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr         0.0009704689
## geom_10m_fl3_L100                                      0.0009548548
## vectorruggedness_hr_ws37                               0.0009332355
## geom_10m_fl4_L140.1                                    0.0008394995
## vectorruggedness_hr_ws57                               0.0008151927
## TRI_hr_ws47                                            0.0008085507
## vectorruggedness_hr_ws35                               0.0007625530
##         
## altpreds  1  2  3  4  5
##        1 26  0  0  0  0
##        2  0 45  0  0  0
##        3  0  0 12  0  0
##        4  0  0  0 18  0
##        5  0  0  0  0  7
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.62962962962963"
## [1] "OOB-error:  0.555555555555556  for predictors geom_10m_fl4_L3 AND TRI_hr_ws35"
##      
## preds  1  2  3  4  5
##     1 15  2  3  2  2
##     2 11 43  9  7  3
##     3  0  0  0  0  0
##     4  0  0  0  9  0
##     5  0  0  0  0  2
## [1] "Kappa overall =  0.444253859348199"
## [1] "Tau overall =  0.548611111111111"
## [1] "mean quality =  0.35752380952381"
## [1] "The quality is  0.35752380952381"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  NaN"
## [1] "classification error rate with altdata:  0.361111111111111"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.601851851851852"
## [1] "OOB-error:  0.601851851851852  for predictors geom_10m_fl4_L3 AND TRI_hr_ws35 AND geom_10m_fl10_L8 AND geom_10m_fl10_L14 AND geom_10m_fl10_L12 AND VerticalDistancetoChannelNetwork_10m AND vectorruggedness_hr_ws55 AND Vertical_Distance_to_Channel_Network_50m AND vectorruggedness_hr_ws27 AND elev_ws15_10m AND geom_10m_fl10_L6 AND DiurnalAnisotropicHeating_10m AND geom_dtm_10m_hyd_fl5_L6 AND Channel_Network_Base_Level_hr AND dtm_hr_CONVEX_r5 AND geom_dtm_10m_hyd_fl5_L3 AND Convexity_10m AND crosc_DTM_50m_avg_ws7_50m AND crosc_ws23_hr_hr AND geom_10m_fl2_L4 AND terraintexture_hr_ws53_tp5 AND slope_ws5_hr_hr AND Plan_Curvature_hr AND minic_ws11_10m AND terraintexture_hr_ws37_tp5 AND terraintexture_hr_ws41_tp5 AND DiurnalAnisotropicHeating_10m AND terraintexture_hr_ws29_tp25 AND Relative_Slope_Position_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr"
##      
## preds  1  2  3  4  5
##     1 26  0  0  0  0
##     2  0 45  0  0  0
##     3  0  0 12  0  0
##     4  0  0  0 18  0
##     5  0  0  0  0  7
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.1 - Minimum precipitation retention capacity"
## [1] "Prediction error at end is:  0.738528138528138"
## [2] "Prediction error at end is:  0.74025974025974" 
## [3] "Prediction error at end is:  0.684848484848485"
## [4] "Prediction error at end is:  0.712554112554112"
## [5] "Prediction error at end is:  0.769264069264069"
## [6] "Prediction error at end is:  0.723376623376623"
##                                               k 1
## 1  geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## 2                                maxic_ws13_hr_hr
## 3 geom_hr_L50m_fl1_rplippatchnum_UE_hr_60cells_hr
## 4                               geom_10m_fl10_L17
## 5                                 geom_10m_fl1_L5
## 6    SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
##                                                 k 2
## 1                       terraintexture_hr_ws53_tp25
## 2      SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 3          Vertical_Distance_to_Channel_Network_50m
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr
## 5     geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr
## 6 geom_hr_L3_fl1_rplippatchdensity_UE_hr_40cells_hr
##                                            k 3
## 1                  terraintexture_hr_ws57_tp25
## 2                  terraintexture_hr_ws53_tp25
## 3                     geom_dtm_10m_hyd_fl5_L30
## 4 geom_hr_L3_fl10_rplipsimpson_UE_hr_5cells_hr
## 5   geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr
## 6                  terraintexture_hr_ws17_tp25
##                                              k 4
## 1                                  profc_ws7_10m
## 2                      vectorstrength_hr_ws17_hr
## 3                    terraintexture_hr_ws53_tp25
## 4 geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr
## 5                                 minic_ws15_10m
## 6                    terraintexture_hr_ws33_tp25
##                                             k 5
## 1                    terraintexture_hr_ws57_tp5
## 2  geom_hr_L3_fl10_rplipshannon_UE_hr_5cells_hr
## 3                              geom_10m_fl10_L9
## 4                    terraintexture_hr_ws41_tp5
## 5                              geom_10m_fl4_L17
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr
## [1] "OBB error with all predictors of  allpreds is  0.675925925925926"
##                             MeanDecreaseGini MeanDecreaseAccuracy
## vectorruggedness_hr_ws11           0.4077027         0.0023059811
## elev_ws5_10m                       0.4122884         0.0011638109
## elev_ws11_10m                      0.3731010         0.0011110040
## vectorruggedness_hr_ws15           0.3706686         0.0010547322
## vectorruggedness_hr_ws23           0.2011323         0.0010062862
## terraintexture_hr_ws57_tp25        0.3201212         0.0009265658
## TotalCurvature_10m                 0.4572167         0.0009247516
## vectorruggedness_hr_ws17           0.2563827         0.0009114961
## Total_Curvature_50m                0.3161605         0.0008903797
## elev_ws3_10m                       0.3370111         0.0008665002
##         
## altpreds  1  2  3  4  5
##        1 12  0  0  0  0
##        2  0 26  0  0  0
##        3  0  0 22  0  0
##        4  0  0  0 26  0
##        5  0  0  0  0 22
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.703703703703704"
## [1] "OOB-error:  0.601851851851852  for predictors geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND maxic_ws13_hr_hr"
##      
## preds  1  2  3  4  5
##     1 10  0  0  0  0
##     2  1 25  0  1  1
##     3  0  0 21  1  0
##     4  0  1  0 23  1
##     5  1  0  1  1 20
## [1] "Kappa overall =  0.894071490845684"
## [1] "Tau overall =  0.895833333333333"
## [1] "mean quality =  0.845974869708003"
## [1] "The quality is  0.845974869708003"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.898014007665387"
## [1] "classification error rate with altdata:  0.0833333333333333"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.685185185185185"
## [1] "OOB-error:  0.537037037037037  for predictors geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND maxic_ws13_hr_hr AND geom_hr_L50m_fl1_rplippatchnum_UE_hr_60cells_hr AND geom_10m_fl10_L17 AND geom_10m_fl1_L5 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND terraintexture_hr_ws53_tp25 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND Vertical_Distance_to_Channel_Network_50m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_40cells_hr AND terraintexture_hr_ws57_tp25 AND terraintexture_hr_ws53_tp25 AND geom_dtm_10m_hyd_fl5_L30 AND geom_hr_L3_fl10_rplipsimpson_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr AND terraintexture_hr_ws17_tp25 AND profc_ws7_10m AND vectorstrength_hr_ws17_hr AND terraintexture_hr_ws53_tp25 AND geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr AND minic_ws15_10m AND terraintexture_hr_ws33_tp25 AND terraintexture_hr_ws57_tp5 AND geom_hr_L3_fl10_rplipshannon_UE_hr_5cells_hr AND geom_10m_fl10_L9 AND terraintexture_hr_ws41_tp5 AND geom_10m_fl4_L17 AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr"
##      
## preds  1  2  3  4  5
##     1 12  0  0  0  0
##     2  0 26  0  0  0
##     3  0  0 22  0  0
##     4  0  0  0 26  0
##     5  0  0  0  0 22
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.2 - Retention capacity for heavy precipitation events"
## [1] "Prediction error at end is:  0.333333333333333"
## [2] "Prediction error at end is:  0.314718614718615"
## [3] "Prediction error at end is:  0.296536796536797"
## [4] "Prediction error at end is:  0.315584415584416"
## [5] "Prediction error at end is:  0.287445887445887"
## [6] "Prediction error at end is:  0.296969696969697"
##                                                k 1
## 1                                 geom_10m_fl8_L16
## 2                      Modified_Catchment_Area_50m
## 3                              Catchment_slope_50m
## 4 geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr
## 5                              Catchment_Area2_50m
## 6                              Total_Curvature_50m
##                                              k 2
## 1                                geom_10m_fl1_L8
## 2                      crosc_DTM_50m_avg_ws5_50m
## 3                                maxic_ws7_hr_hr
## 4                                longc_ws3_hr_hr
## 5       geom_hr_L3_fl1_rplipmps_UE_hr_60cells_hr
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr
##                                                k 3
## 1         geom_hr_L3_fl10_rplipmps_UE_hr_5cells_hr
## 2                                  geom_10m_fl2_L8
## 3                                   crosc_ws11_10m
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr
## 5                              Catchment_slope_50m
## 6                                  SlopeHeight_10m
##                                            k 4                       k 5
## 1                            geom_10m_fl10_L16          geom_10m_fl3_L11
## 2 geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr geom_dtm_10m_hyd_fl5_L130
## 3     Vertical_Distance_to_Channel_Network_50m         geom_10m_fl3_L120
## 4                             fischerk_ws53_hr           geom_10m_fl1_L9
## 5                        FlowLineCurvature_10m               Texture_50m
## 6                     vectorruggedness_hr_ws35 geom_dtm_10m_hyd_fl5_L110
## [1] "OBB error with all predictors of  allpreds is  0.287037037037037"
##                                          MeanDecreaseGini
## FlowLineCurvature_10m                           0.3883731
## vectorruggedness_hr_ws25                        0.1386517
## crosc_ws11_10m                                  0.2424327
## vectorruggedness_hr_ws21                        0.2543887
## minic_ws29_hr_hr                                0.3022987
## profc_ws15_hr_hr                                0.1493986
## maxic_ws15_hr_hr                                0.1506604
## vectorruggedness_hr_ws47                        0.2067361
## Vertical_Distance_to_Channel_Network_50m        0.2598812
## profc_ws19_hr_hr                                0.1860083
##                                          MeanDecreaseAccuracy
## FlowLineCurvature_10m                             0.003774068
## vectorruggedness_hr_ws25                          0.002113831
## crosc_ws11_10m                                    0.002012501
## vectorruggedness_hr_ws21                          0.001863489
## minic_ws29_hr_hr                                  0.001776280
## profc_ws15_hr_hr                                  0.001477575
## maxic_ws15_hr_hr                                  0.001384619
## vectorruggedness_hr_ws47                          0.001379671
## Vertical_Distance_to_Channel_Network_50m          0.001352361
## profc_ws19_hr_hr                                  0.001246409
##         
## altpreds  1  2  3  4  5
##        1 75  0  0  0  0
##        2  0 11  0  0  0
##        3  0  0  6  0  0
##        4  0  0  0  8  0
##        5  0  0  0  0  8
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.296296296296296"
## [1] "OOB-error:  0.305555555555556  for predictors geom_10m_fl8_L16 AND Modified_Catchment_Area_50m"
##      
## preds  1  2  3  4  5
##     1 75  5  6  2  5
##     2  0  6  0  0  0
##     3  0  0  0  0  0
##     4  0  0  0  6  0
##     5  0  0  0  0  3
## [1] "Kappa overall =  0.572841133816744"
## [1] "Tau overall =  0.791666666666667"
## [1] "mean quality =  0.495381231671554"
## [1] "The quality is  0.495381231671554"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  NaN"
## [1] "classification error rate with altdata:  0.166666666666667"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.296296296296296"
## [1] "OOB-error:  0.25  for predictors geom_10m_fl8_L16 AND Modified_Catchment_Area_50m AND Catchment_slope_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND Catchment_Area2_50m AND Total_Curvature_50m AND geom_10m_fl1_L8 AND crosc_DTM_50m_avg_ws5_50m AND maxic_ws7_hr_hr AND longc_ws3_hr_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipmps_UE_hr_5cells_hr AND geom_10m_fl2_L8 AND crosc_ws11_10m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND Catchment_slope_50m AND SlopeHeight_10m AND geom_10m_fl10_L16 AND geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr AND Vertical_Distance_to_Channel_Network_50m AND fischerk_ws53_hr AND FlowLineCurvature_10m AND vectorruggedness_hr_ws35 AND geom_10m_fl3_L11 AND geom_dtm_10m_hyd_fl5_L130 AND geom_10m_fl3_L120 AND geom_10m_fl1_L9 AND Texture_50m AND geom_dtm_10m_hyd_fl5_L110"
##      
## preds  1  2  3  4  5
##     1 75  0  0  0  0
##     2  0 11  0  0  0
##     3  0  0  6  0  0
##     4  0  0  0  8  0
##     5  0  0  0  0  8
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.3 - groundwater reformation rate"
## [1] "Prediction error at end is:  0.704761904761905"
## [2] "Prediction error at end is:  0.731601731601732"
## [3] "Prediction error at end is:  0.665367965367965"
## [4] "Prediction error at end is:  0.654978354978355"
## [5] "Prediction error at end is:  0.598701298701299"
## [6] "Prediction error at end is:  0.637229437229437"
##                                               k 1
## 1                                 geom_10m_fl8_L5
## 2  geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr
## 3 geom_hr_L50m_fl1_rpliprichness_UE_hr_10cells_hr
## 4                                 geom_10m_fl8_L4
## 5                                 geom_10m_fl8_L3
## 6  geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr
##                                            k 2
## 1                                          SGU
## 2 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 3                    slope_DTM_50m_avg_ws7_50m
## 4                                  SGUT_wTGnew
## 5                    maxic_DTM_50m_avg_ws5_50m
## 6                             geom_10m_fl10_L6
##                                        k 3
## 1                                      SGU
## 2                         geom_10m_fl8_L50
## 3               CrossSectionalCurvature_hr
## 4 Vertical_Distance_to_Channel_Network_50m
## 5               terraintexture_hr_ws29_tp5
## 6                       TotalCurvature_10m
##                                         k 4
## 1               terraintexture_hr_ws49_tp25
## 2                               TRI_hr_ws17
## 3 geom_hr_L3_fl10_rplipmps_UE_hr_20cells_hr
## 4                          fischerk_ws33_hr
## 5                             profc_ws7_10m
## 6                terraintexture_hr_ws49_tp5
##                                                k 5
## 1                      terraintexture_hr_ws37_tp25
## 2 geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr
## 3                                  geom_10m_fl8_L7
## 4                        vectorstrength_hr_ws13_hr
## 5                                 geom_10m_fl2_L16
## 6                                geom_10m_fl1_L110
## [1] "OBB error with all predictors of  allpreds is  0.62962962962963"
##                          MeanDecreaseGini MeanDecreaseAccuracy
## vectorruggedness_hr_ws59        0.3488135          0.001858079
## vectorruggedness_hr_ws43        0.3607745          0.001646025
## vectorruggedness_hr_ws35        0.3500364          0.001586086
## vectorruggedness_hr_ws55        0.4023994          0.001418956
## vectorruggedness_hr_ws13        0.3501757          0.001317820
## vectorruggedness_hr_ws27        0.3841066          0.001313411
## LSFactor_10m                    0.1499655          0.001205072
## vectorruggedness_hr_ws51        0.2968573          0.001154607
## SGU.1                           0.4128421          0.001094913
## TRI_hr_ws5                      0.2039083          0.001091161
##         
## altpreds  1  2  3  4  5
##        1 14  0  0  0  0
##        2  0 39  0  0  0
##        3  0  0 14  0  0
##        4  0  0  0 19  0
##        5  0  0  0  0 22
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.62037037037037"
## [1] "OOB-error:  0.509259259259259  for predictors geom_10m_fl8_L5 AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr"
##      
## preds  1  2  3  4  5
##     1  6  0  3  0  1
##     2  4 34  5  2  5
##     3  0  0  5  0  0
##     4  0  0  0 15  0
##     5  4  5  1  2 16
## [1] "Kappa overall =  0.598279669882599"
## [1] "Tau overall =  0.62962962962963"
## [1] "mean quality =  0.513743985632531"
## [1] "The quality is  0.513743985632531"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.640367987990471"
## [1] "classification error rate with altdata:  0.296296296296296"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.675925925925926"
## [1] "OOB-error:  0.546296296296296  for predictors geom_10m_fl8_L5 AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl8_L4 AND geom_10m_fl8_L3 AND geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr AND SGU AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND slope_DTM_50m_avg_ws7_50m AND SGUT_wTGnew AND maxic_DTM_50m_avg_ws5_50m AND geom_10m_fl10_L6 AND SGU AND geom_10m_fl8_L50 AND CrossSectionalCurvature_hr AND Vertical_Distance_to_Channel_Network_50m AND terraintexture_hr_ws29_tp5 AND TotalCurvature_10m AND terraintexture_hr_ws49_tp25 AND TRI_hr_ws17 AND geom_hr_L3_fl10_rplipmps_UE_hr_20cells_hr AND fischerk_ws33_hr AND profc_ws7_10m AND terraintexture_hr_ws49_tp5 AND terraintexture_hr_ws37_tp25 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl8_L7 AND vectorstrength_hr_ws13_hr AND geom_10m_fl2_L16 AND geom_10m_fl1_L110"
##      
## preds  1  2  3  4  5
##     1 14  0  0  0  0
##     2  0 39  0  0  0
##     3  0  0 14  0  0
##     4  0  0  0 19  0
##     5  0  0  0  0 22
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.4 - Potential for providing nutrients for plants"
## [1] "Prediction error at end is:  0.25974025974026" 
## [2] "Prediction error at end is:  0.333766233766234"
## [3] "Prediction error at end is:  0.315584415584416"
## [4] "Prediction error at end is:  0.296536796536797"
## [5] "Prediction error at end is:  0.306493506493506"
## [6] "Prediction error at end is:  0.324675324675325"
##                                                  k 1
## 1                                                SGU
## 2                                   geom_10m_fl10_L3
## 3                          vectorstrength_hr_ws17_hr
## 4                                    geom_10m_fl2_L5
## 5 geom_hr_L50m_fl1_rplipedgedensity_UE_hr_60cells_hr
## 6               VerticalDistancetoChannelNetwork_10m
##                                                  k 2
## 1                                                SGU
## 2                                   geom_10m_fl8_L11
## 3 geom_hr_L3_fl10_rplippatchdensity_UE_hr_60cells_hr
## 4        geom_hr_L3_fl1_rplipsimpson_UE_hr_5cells_hr
## 5                                    geom_10m_fl2_L4
## 6                  sagaTopographic_Wetness_Index_50m
##                                             k 3
## 1                                           SGU
## 2                    Longitudinal_Curvature_50m
## 3                               geom_10m_fl1_L9
## 4 geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr
## 5      geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr
## 6                               geom_10m_fl1_L4
##                                                   k 4
## 1                                                 SGU
## 2                                     ValleyDepth_10m
## 3                VerticalDistancetoChannelNetwork_10m
## 4                                         Texture_10m
## 5                         terraintexture_hr_ws29_tp25
## 6 geom_hr_L50m_fl1_rplippatchdensity_UE_hr_10cells_hr
##                                       k 5
## 1                                     SGU
## 2                        geom_10m_fl3_L12
## 3 Vertical_Distance_to_Channel_Network_hr
## 4                     Catchment_slope_50m
## 5                        geom_10m_fl1_L30
## 6                        geom_10m_fl3_L23
## [1] "OBB error with all predictors of  allpreds is  0.314814814814815"
##                          MeanDecreaseGini MeanDecreaseAccuracy
## SGU                             0.6680496          0.004606939
## vectorruggedness_hr_ws25        0.6748379          0.004380943
## SGU.1                           0.6701422          0.004022976
## vectorruggedness_hr_ws27        0.4492908          0.003844611
## minic_ws7_hr_hr                 0.3908256          0.002601657
## TotalCurvature_10m              0.5374482          0.002522411
## vectorruggedness_hr_ws29        0.5675637          0.002382132
## vectorruggedness_hr_ws31        0.3422939          0.002213987
## elev_ws11_10m                   0.3290646          0.002190745
## vectorruggedness_hr_ws13        0.2573725          0.002075946
##         
## altpreds  1  3  5
##        1 65  0  0
##        3  0 11  0
##        5  0  0 32
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.305555555555556"
## [1] "OOB-error:  0.305555555555556  for predictors SGU AND geom_10m_fl10_L3"
##      
## preds  1  3  5
##     1 60  4  4
##     3  0  1  0
##     5  5  6 28
## [1] "Kappa overall =  0.657142857142857"
## [1] "Tau overall =  0.736111111111111"
## [1] "mean quality =  0.521329896608648"
## [1] "The quality is  0.521329896608648"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.572763122715766"
## [1] "classification error rate with altdata:  0.175925925925926"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.296296296296296"
## [1] "OOB-error:  0.231481481481481  for predictors SGU AND geom_10m_fl10_L3 AND vectorstrength_hr_ws17_hr AND geom_10m_fl2_L5 AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_60cells_hr AND VerticalDistancetoChannelNetwork_10m AND SGU AND geom_10m_fl8_L11 AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_5cells_hr AND geom_10m_fl2_L4 AND sagaTopographic_Wetness_Index_50m AND SGU AND Longitudinal_Curvature_50m AND geom_10m_fl1_L9 AND geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND geom_10m_fl1_L4 AND SGU AND ValleyDepth_10m AND VerticalDistancetoChannelNetwork_10m AND Texture_10m AND terraintexture_hr_ws29_tp25 AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_10cells_hr AND SGU AND geom_10m_fl3_L12 AND Vertical_Distance_to_Channel_Network_hr AND Catchment_slope_50m AND geom_10m_fl1_L30 AND geom_10m_fl3_L23"
##      
## preds  1  3  5
##     1 65  0  0
##     3  0 11  0
##     5  0  0 32
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.5 - Potential as a CO2 sink"
## [1] "Prediction error at end is:  0.541991341991342"
## [2] "Prediction error at end is:  0.466666666666667"
## [3] "Prediction error at end is:  0.42987012987013" 
## [4] "Prediction error at end is:  0.382251082251082"
## [5] "Prediction error at end is:  0.410822510822511"
## [6] "Prediction error at end is:  0.392640692640693"
##                                            k 1
## 1                             slope_ws19_hr_hr
## 2                             fischerk_ws19_hr
## 3                    vectorstrength_hr_ws27_hr
## 4 geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr
## 5                   minic_DTM_50m_avg_ws11_50m
## 6                                  TRI_hr_ws18
##                                              k 2
## 1                               slope_ws19_hr_hr
## 2                      maxic_DTM_50m_avg_ws5_50m
## 3                        vectorruggedness_hr_ws7
## 4     geom_hr_L3_fl1_rplipshape_UE_hr_40cells_hr
## 5                              geom_10m_fl3_L110
## 6 geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr
##                                                k 3
## 1                       CrossSectionalCurvature_hr
## 2                               TotalCurvature_10m
## 3                    Topographic_Wetness_Index_50m
## 4                    Channel_Network_Base_Level_hr
## 5 geom_hr_L50m_fl1_rplipdominance_UE_hr_20cells_hr
## 6                                  longc_ws5_hr_hr
##                                            k 4                      k 5
## 1                              crosc_ws5_hr_hr                      SGU
## 2                             slope_ws23_hr_hr vectorruggedness_hr_ws35
## 3                     geom_dtm_10m_hyd_fl5_L29     MassBalanceIndex_10m
## 4 geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr            crosc_ws3_10m
## 5                   CrossSectionalCurvature_hr              SGUT_wTGnew
## 6                            geom_10m_fl2_L140         geom_10m_fl4_L19
## [1] "OBB error with all predictors of  allpreds is  0.333333333333333"
##                          MeanDecreaseGini MeanDecreaseAccuracy
## elev_ws3_10m                    0.9418425          0.007665726
## vectorruggedness_hr_ws7         0.7755776          0.006588824
## elev_ws15_10m                   0.7687863          0.006467315
## vectorruggedness_hr_ws9         0.3699884          0.005604114
## vectorruggedness_hr_ws21        0.5762638          0.005282938
## vectorruggedness_hr_ws3         0.8258676          0.005266290
## elev_ws11_10m                   0.6622378          0.005152287
## elev_ws5_10m                    0.6380886          0.005086209
## vectorruggedness_hr_ws5         0.6236385          0.004842494
## vectorruggedness_hr_ws31        0.7160756          0.004086461
##         
## altpreds  1  2  3  4  5
##        1 51  0  0  0  0
##        2  0  1  0  0  0
##        3  0  0 18  0  0
##        4  0  0  0 33  0
##        5  0  0  0  0  5
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.305555555555556"
## [1] "OOB-error:  0.314814814814815  for predictors slope_ws19_hr_hr AND fischerk_ws19_hr"
##      
## preds  1  2  3  4  5
##     1 51  0  0  0  0
##     2  0  1  0  0  0
##     3  0  0 18  0  0
##     4  0  0  0 33  0
##     5  0  0  0  0  5
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.333333333333333"
## [1] "OOB-error:  0.240740740740741  for predictors slope_ws19_hr_hr AND fischerk_ws19_hr AND vectorstrength_hr_ws27_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr AND minic_DTM_50m_avg_ws11_50m AND TRI_hr_ws18 AND slope_ws19_hr_hr AND maxic_DTM_50m_avg_ws5_50m AND vectorruggedness_hr_ws7 AND geom_hr_L3_fl1_rplipshape_UE_hr_40cells_hr AND geom_10m_fl3_L110 AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND CrossSectionalCurvature_hr AND TotalCurvature_10m AND Topographic_Wetness_Index_50m AND Channel_Network_Base_Level_hr AND geom_hr_L50m_fl1_rplipdominance_UE_hr_20cells_hr AND longc_ws5_hr_hr AND crosc_ws5_hr_hr AND slope_ws23_hr_hr AND geom_dtm_10m_hyd_fl5_L29 AND geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr AND CrossSectionalCurvature_hr AND geom_10m_fl2_L140 AND SGU AND vectorruggedness_hr_ws35 AND MassBalanceIndex_10m AND crosc_ws3_10m AND SGUT_wTGnew AND geom_10m_fl4_L19"
##      
## preds  1  2  3  4  5
##     1 51  0  0  0  0
##     2  0  1  0  0  0
##     3  0  0 18  0  0
##     4  0  0  0 33  0
##     5  0  0  0  0  5
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.1 - Potential for retention of heavy metals"
## [1] "Prediction error at end is:  0.444588744588745"
## [2] "Prediction error at end is:  0.481385281385281"
## [3] "Prediction error at end is:  0.407792207792208"
## [4] "Prediction error at end is:  0.398268398268398"
## [5] "Prediction error at end is:  0.361471861471861"
## [6] "Prediction error at end is:  0.38961038961039" 
##                             k 1
## 1                           SGU
## 2      vectorruggedness_hr_ws45
## 3 Channel_Network_Base_Level_hr
## 4           Catchment_slope_50m
## 5              geom_10m_fl4_L47
## 6              fischerk_ws17_hr
##                                                k 2
## 1                                 geom_10m_fl8_L60
## 2                                    Convexity_50m
## 3                                 geom_10m_fl8_L50
## 4                      Relative_Slope_Position_50m
## 5 geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr
## 6                                 geom_10m_fl4_L14
##                                            k 3
## 1                              geom_10m_fl8_L8
## 2 geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr
## 3                     geom_dtm_10m_hyd_fl5_L60
## 4                             geom_10m_fl4_L23
## 5 geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr
## 6                              Valley_Depth_hr
##                                              k 4
## 1                               geom_10m_fl8_L12
## 2 geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr
## 3                               geom_10m_fl3_L25
## 4                      profc_DTM_50m_avg_ws7_50m
## 5                               geom_10m_fl1_L50
## 6  geom_hr_L3_fl10_rplipshannon_UE_hr_40cells_hr
##                                        k 5
## 1                                      SGU
## 2                         geom_10m_fl1_L42
## 3                  TangentialCurvature_10m
## 4              ChannelNetworkBaseLevel_10m
## 5 geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr
## 6                 vectorruggedness_hr_ws59
## [1] "OBB error with all predictors of  allpreds is  0.342592592592593"
##                                          MeanDecreaseGini
## vectorruggedness_hr_ws7                         0.3681345
## elev_ws5_10m                                    0.4273979
## vectorruggedness_hr_ws23                        0.3559804
## vectorruggedness_hr_ws41                        0.5187288
## elev_ws3_10m                                    0.7682106
## vectorruggedness_hr_ws39                        0.5674367
## vectorruggedness_hr_ws29                        0.5086838
## Vertical_Distance_to_Channel_Network_50m        0.4177763
## vectorruggedness_hr_ws55                        0.4823247
## vectorruggedness_hr_ws25                        0.5943780
##                                          MeanDecreaseAccuracy
## vectorruggedness_hr_ws7                           0.003925552
## elev_ws5_10m                                      0.003555039
## vectorruggedness_hr_ws23                          0.003483173
## vectorruggedness_hr_ws41                          0.003437154
## elev_ws3_10m                                      0.003407096
## vectorruggedness_hr_ws39                          0.003403081
## vectorruggedness_hr_ws29                          0.003312758
## Vertical_Distance_to_Channel_Network_50m          0.003246063
## vectorruggedness_hr_ws55                          0.003198155
## vectorruggedness_hr_ws25                          0.003189255
##         
## altpreds  1  2  3  4  5
##        1 42  0  0  0  0
##        2  0  7  0  0  0
##        3  0  0 10  0  0
##        4  0  0  0  7  0
##        5  0  0  0  0 42
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.333333333333333"
## [1] "OOB-error:  0.333333333333333  for predictors SGU AND vectorruggedness_hr_ws45"
##      
## preds  1  2  3  4  5
##     1 41  2  3  2  5
##     2  0  5  0  0  0
##     3  0  0  6  0  0
##     4  0  0  0  3  0
##     5  1  0  1  2 37
## [1] "Kappa overall =  0.772631578947368"
## [1] "Tau overall =  0.814814814814815"
## [1] "mean quality =  0.661292845640672"
## [1] "The quality is  0.661292845640672"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.768075773614234"
## [1] "classification error rate with altdata:  0.148148148148148"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.342592592592593"
## [1] "OOB-error:  0.314814814814815  for predictors SGU AND vectorruggedness_hr_ws45 AND Channel_Network_Base_Level_hr AND Catchment_slope_50m AND geom_10m_fl4_L47 AND fischerk_ws17_hr AND geom_10m_fl8_L60 AND Convexity_50m AND geom_10m_fl8_L50 AND Relative_Slope_Position_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr AND geom_10m_fl4_L14 AND geom_10m_fl8_L8 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr AND geom_dtm_10m_hyd_fl5_L60 AND geom_10m_fl4_L23 AND geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr AND Valley_Depth_hr AND geom_10m_fl8_L12 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr AND geom_10m_fl3_L25 AND profc_DTM_50m_avg_ws7_50m AND geom_10m_fl1_L50 AND geom_hr_L3_fl10_rplipshannon_UE_hr_40cells_hr AND SGU AND geom_10m_fl1_L42 AND TangentialCurvature_10m AND ChannelNetworkBaseLevel_10m AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND vectorruggedness_hr_ws59"
##      
## preds  1  2  3  4  5
##     1 42  0  0  0  0
##     2  0  7  0  0  0
##     3  0  0 10  0  0
##     4  0  0  0  7  0
##     5  0  0  0  0 42
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.2 - Potential for transforming organic contaminants"
## [1] "Prediction error at end is:  0.6"              
## [2] "Prediction error at end is:  0.572294372294372"
## [3] "Prediction error at end is:  0.654978354978355"
## [4] "Prediction error at end is:  0.617748917748918"
## [5] "Prediction error at end is:  0.627705627705628"
## [6] "Prediction error at end is:  0.637662337662338"
##                                              k 1
## 1                              geom_10m_fl10_L21
## 2                      geom_dtm_10m_hyd_fl5_L110
## 3 geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr
## 4                     minic_DTM_50m_avg_ws11_50m
## 5                             geom_10m_fl10_L150
## 6                     maxic_DTM_50m_avg_ws11_50m
##                                              k 2                       k 3
## 1                        geom_dtm_10m_hyd_fl5_L8         geom_10m_fl10_L90
## 2   geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr              elev_ws3_10m
## 3  geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr planc_DTM_50m_avg_ws3_50m
## 4 geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr        geom_10m_fl10_L120
## 5  geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr           maxic_ws3_hr_hr
## 6  geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr         geom_10m_fl8_L120
##                                              k 4
## 1                              geom_10m_fl10_L22
## 2                              PlanCurvature_10m
## 3 geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr
## 4                                geom_10m_fl3_L6
## 5                              geom_10m_fl10_L17
## 6                              geom_10m_fl10_L19
##                                              k 5
## 1                              geom_10m_fl10_L22
## 2                                  elev_ws15_10m
## 3                                geom_10m_fl8_L9
## 4                                  Convexity_50m
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr
## 6                                  minic_ws3_10m
## [1] "OBB error with all predictors of  allpreds is  0.527777777777778"
##                                          MeanDecreaseGini
## elev_ws5_10m                                    0.5027974
## Vertical_Distance_to_Channel_Network_50m        0.5260216
## vectorruggedness_hr_ws27                        0.6309566
## vectorruggedness_hr_ws25                        0.4315166
## vectorruggedness_hr_ws33                        0.3605873
## vectorruggedness_hr_ws37                        0.3469170
## elev_ws15_10m                                   0.3954492
## vectorruggedness_hr_ws31                        0.3766044
## Channel_Network_Base_Level_hr                   0.2919898
## elev_ws11_10m                                   0.4945420
##                                          MeanDecreaseAccuracy
## elev_ws5_10m                                      0.003366026
## Vertical_Distance_to_Channel_Network_50m          0.003336541
## vectorruggedness_hr_ws27                          0.002601471
## vectorruggedness_hr_ws25                          0.002287444
## vectorruggedness_hr_ws33                          0.001991645
## vectorruggedness_hr_ws37                          0.001954972
## elev_ws15_10m                                     0.001914333
## vectorruggedness_hr_ws31                          0.001840628
## Channel_Network_Base_Level_hr                     0.001834953
## elev_ws11_10m                                     0.001675054
##         
## altpreds  1  2  3  4  5
##        1 10  0  0  0  0
##        2  0 22  0  0  0
##        3  0  0 28  0  0
##        4  0  0  0  9  0
##        5  0  0  0  0 39
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.537037037037037"
## [1] "OOB-error:  0.5  for predictors geom_10m_fl10_L21 AND geom_dtm_10m_hyd_fl5_L110"
##      
## preds  1  2  3  4  5
##     1  4  1  0  0  1
##     2  4 15  7  2  4
##     3  2  3 16  0  3
##     4  0  0  0  2  0
##     5  0  3  5  5 31
## [1] "Kappa overall =  0.491405698139864"
## [1] "Tau overall =  0.537037037037037"
## [1] "mean quality =  0.396153846153846"
## [1] "The quality is  0.396153846153846"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.507139419578754"
## [1] "classification error rate with altdata:  0.37037037037037"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.555555555555556"
## [1] "OOB-error:  0.453703703703704  for predictors geom_10m_fl10_L21 AND geom_dtm_10m_hyd_fl5_L110 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr AND minic_DTM_50m_avg_ws11_50m AND geom_10m_fl10_L150 AND maxic_DTM_50m_avg_ws11_50m AND geom_dtm_10m_hyd_fl5_L8 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl10_L90 AND elev_ws3_10m AND planc_DTM_50m_avg_ws3_50m AND geom_10m_fl10_L120 AND maxic_ws3_hr_hr AND geom_10m_fl8_L120 AND geom_10m_fl10_L22 AND PlanCurvature_10m AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr AND geom_10m_fl3_L6 AND geom_10m_fl10_L17 AND geom_10m_fl10_L19 AND geom_10m_fl10_L22 AND elev_ws15_10m AND geom_10m_fl8_L9 AND Convexity_50m AND geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr AND minic_ws3_10m"
##      
## preds  1  2  3  4  5
##     1 10  0  0  0  0
##     2  0 22  0  0  0
##     3  0  0 28  0  0
##     4  0  0  0  9  0
##     5  0  0  0  0 39
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"
## [1] "Prediction error at end is:  0.0826839826839827"
## [2] "Prediction error at end is:  0.0463203463203463"
## [3] "Prediction error at end is:  0.0463203463203463"
## [4] "Prediction error at end is:  0.0372294372294372"
## [5] "Prediction error at end is:  0.0372294372294372"
## [6] "Prediction error at end is:  0.0372294372294372"
##                                                 k 1
## 1                                  Slope_Height_50m
## 2                                  geom_10m_fl2_L15
## 3                                  geom_10m_fl2_L16
## 4                                  geom_10m_fl2_L17
## 5                                  geom_10m_fl2_L18
## 6 geom_hr_L3_fl10_rplipedgedensity_UE_hr_40cells_hr
##                                               k 2
## 1 geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr
## 2                       planc_DTM_50m_avg_ws5_50m
## 3                                geom_10m_fl2_L19
## 4                                geom_10m_fl2_L22
## 5                                 geom_10m_fl4_L3
## 6                                geom_10m_fl8_L17
##                                               k 3
## 1                                planc_ws19_hr_hr
## 2                            MinimalCurvature_10m
## 3                                crosc_ws11_hr_hr
## 4 geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr
## 5 geom_hr_L50m_fl1_rpliprichness_UE_hr_60cells_hr
## 6                                fischerk_ws31_hr
##                                                k 4
## 1                                  maxic_ws7_hr_hr
## 2 geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr
## 3 geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr
## 4                                 geom_10m_fl10_L3
## 5                                 minic_ws15_hr_hr
## 6                                geom_10m_fl10_L11
##                                               k 5
## 1 geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr
## 2                                geom_10m_fl2_L19
## 3                                 geom_10m_fl4_L8
## 4                                geom_10m_fl1_L60
## 5                                geom_10m_fl2_L60
## 6                                   minic_ws7_10m
## [1] "OBB error with all predictors of  allpreds is  0.037037037037037"
##                                                   MeanDecreaseGini
## geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr         0.046223453
## geom_10m_fl2_L17.1                                     0.082940315
## geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr           0.029509434
## PlanCurvature_10m                                      0.054469869
## Mid_Slope_Positon_50m                                  0.085348233
## geom_10m_fl3_L130.1                                    0.035999554
## geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr             0.071236465
## geom_10m_fl8_L40                                       0.003010885
## dtm_hr_TEXTURE_r20                                     0.007466083
## geom_hr_L3_fl1_rplippatchdensity_UE_hr_10cells_hr      0.006166667
##                                                   MeanDecreaseAccuracy
## geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr            0.0004327331
## geom_10m_fl2_L17.1                                        0.0004074340
## geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr              0.0003541687
## PlanCurvature_10m                                         0.0003217833
## Mid_Slope_Positon_50m                                     0.0003035496
## geom_10m_fl3_L130.1                                       0.0002844969
## geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr                0.0002628852
## geom_10m_fl8_L40                                          0.0002554463
## dtm_hr_TEXTURE_r20                                        0.0002481203
## geom_hr_L3_fl1_rplippatchdensity_UE_hr_10cells_hr         0.0002417582
##         
## altpreds   4   5
##        4   4   0
##        5   0 104
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.037037037037037"
## [1] "OOB-error:  0.0185185185185185  for predictors Slope_Height_50m AND geom_10m_fl2_L15"
##      
## preds   4   5
##     4   2   0
##     5   2 104
## [1] "Kappa overall =  0.658227848101266"
## [1] "Tau overall =  0.962962962962963"
## [1] "mean quality =  0.740566037735849"
## [1] "The quality is  0.740566037735849"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.518568491248852"
## [1] "classification error rate with altdata:  0.0185185185185185"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.037037037037037"
## [1] "OOB-error:  0.0277777777777778  for predictors Slope_Height_50m AND geom_10m_fl2_L15 AND geom_10m_fl2_L16 AND geom_10m_fl2_L17 AND geom_10m_fl2_L18 AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_40cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND planc_DTM_50m_avg_ws5_50m AND geom_10m_fl2_L19 AND geom_10m_fl2_L22 AND geom_10m_fl4_L3 AND geom_10m_fl8_L17 AND planc_ws19_hr_hr AND MinimalCurvature_10m AND crosc_ws11_hr_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_60cells_hr AND fischerk_ws31_hr AND maxic_ws7_hr_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND geom_10m_fl10_L3 AND minic_ws15_hr_hr AND geom_10m_fl10_L11 AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND geom_10m_fl2_L19 AND geom_10m_fl4_L8 AND geom_10m_fl1_L60 AND geom_10m_fl2_L60 AND minic_ws7_10m"
##      
## preds   4   5
##     4   4   0
##     5   0 104
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  0.870192307692308"
## [1] "classification error rate with altdata:  0"
## [1] "1d.4 - Potential for retention of water-soluble contaminants"
## [1] "Prediction error at end is:  0.583549783549784"
## [2] "Prediction error at end is:  0.582683982683983"
## [3] "Prediction error at end is:  0.601298701298701"
## [4] "Prediction error at end is:  0.583116883116883"
## [5] "Prediction error at end is:  0.582683982683983"
## [6] "Prediction error at end is:  0.564069264069264"
##                           k 1
## 1             geom_10m_fl8_L5
## 2 Relative_Slope_Position_50m
## 3   minic_DTM_50m_avg_ws5_50m
## 4            geom_10m_fl2_L15
## 5             geom_10m_fl1_L7
## 6           Plan_Curvature_hr
##                                                 k 2
## 1                                 geom_10m_fl3_L120
## 2    geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr
## 3    geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr
## 4                                 geom_10m_fl3_L130
## 5        geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr
## 6 geom_hr_L50m_fl1_rplipedgedensity_UE_hr_5cells_hr
##                           k 3                                          k 4
## 1             geom_10m_fl8_L4                            geom_10m_fl8_L110
## 2 Relative_Slope_Position_50m                   longc_DTM_50m_avg_ws11_50m
## 3                 TRI_hr_ws53                     vectorruggedness_hr_ws35
## 4             maxic_ws5_hr_hr                               slope_ws15_10m
## 5             geom_10m_fl8_L5 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 6        NormalizedHeight_10m                              profc_ws5_hr_hr
##                                                 k 5
## 1                                 geom_10m_fl10_L14
## 2                                 geom_10m_fl10_L10
## 3     geom_hr_L3_fl1_rpliprichness_UE_hr_40cells_hr
## 4  geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr
## 5                          geom_dtm_10m_hyd_fl5_L32
## 6 geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr
## [1] "OBB error with all predictors of  allpreds is  0.5"
##                                   MeanDecreaseGini MeanDecreaseAccuracy
## sagaTopographic_Wetness_Index_50m        0.4100743          0.002946161
## vectorruggedness_hr_ws35                 0.4577494          0.002292830
## vectorruggedness_hr_ws41                 0.4793662          0.002267833
## SGU.1                                    0.2898393          0.002259740
## TRI_hr_ws28                              0.2288264          0.001755498
## elev_ws11_10m                            0.4273840          0.001746943
## vectorruggedness_hr_ws25                 0.3132226          0.001589001
## vectorruggedness_hr_ws23                 0.3064092          0.001528978
## vectorruggedness_hr_ws27                 0.3221944          0.001405396
## TWI_10m                                  0.1942418          0.001397361
##         
## altpreds  1  2  3  4  5
##        1 47  0  0  0  0
##        2  0 14  0  0  0
##        3  0  0 17  0  0
##        4  0  0  0 22  0
##        5  0  0  0  0  8
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.481481481481481"
## [1] "OOB-error:  0.425925925925926  for predictors geom_10m_fl8_L5 AND Relative_Slope_Position_50m"
##      
## preds  1  2  3  4  5
##     1 46  4  4  4  0
##     2  0  7  1  1  0
##     3  0  0  7  0  0
##     4  1  3  5 17  3
##     5  0  0  0  0  5
## [1] "Kappa overall =  0.649656893325016"
## [1] "Tau overall =  0.699074074074074"
## [1] "mean quality =  0.550785144566301"
## [1] "The quality is  0.550785144566301"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.675892376003383"
## [1] "classification error rate with altdata:  0.240740740740741"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.518518518518518"
## [1] "OOB-error:  0.5  for predictors geom_10m_fl8_L5 AND Relative_Slope_Position_50m AND minic_DTM_50m_avg_ws5_50m AND geom_10m_fl2_L15 AND geom_10m_fl1_L7 AND Plan_Curvature_hr AND geom_10m_fl3_L120 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr AND geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr AND geom_10m_fl3_L130 AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl8_L4 AND Relative_Slope_Position_50m AND TRI_hr_ws53 AND maxic_ws5_hr_hr AND geom_10m_fl8_L5 AND NormalizedHeight_10m AND geom_10m_fl8_L110 AND longc_DTM_50m_avg_ws11_50m AND vectorruggedness_hr_ws35 AND slope_ws15_10m AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND profc_ws5_hr_hr AND geom_10m_fl10_L14 AND geom_10m_fl10_L10 AND geom_hr_L3_fl1_rpliprichness_UE_hr_40cells_hr AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND geom_dtm_10m_hyd_fl5_L32 AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr"
##      
## preds  1  2  3  4  5
##     1 47  0  0  0  0
##     2  0 14  0  0  0
##     3  0  0 17  0  0
##     4  0  0  0 22  0
##     5  0  0  0  0  8
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.5 - Potential as buffer for acidic contaminants"
## [1] "Prediction error at end is:  0.676190476190476"
## [2] "Prediction error at end is:  0.722077922077922"
## [3] "Prediction error at end is:  0.712987012987013"
## [4] "Prediction error at end is:  0.703030303030303"
## [5] "Prediction error at end is:  0.667099567099567"
## [6] "Prediction error at end is:  0.603030303030303"
##                           k 1                       k 2
## 1   geom_dtm_10m_hyd_fl5_L110 geom_dtm_10m_hyd_fl5_L140
## 2          Catchment_area_50m planc_DTM_50m_avg_ws9_50m
## 3            geom_10m_fl4_L40                       SGU
## 4             planc_ws5_hr_hr           crosc_ws5_hr_hr
## 5    geom_dtm_10m_hyd_fl5_L42          geom_10m_fl10_L3
## 6 terraintexture_hr_ws53_tp25   geom_dtm_10m_hyd_fl5_L4
##                                             k 3
## 1 geom_hr_L3_fl10_rplipshannon_UE_hr_10cells_hr
## 2                     geom_dtm_10m_hyd_fl5_L110
## 3                               maxic_ws3_hr_hr
## 4                    terraintexture_hr_ws41_tp5
## 5                              planc_ws29_hr_hr
## 6                                 LS_Factor_50m
##                                                  k 4                   k 5
## 1                                   geom_10m_fl1_L12       geom_10m_fl1_L8
## 2 geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr      geom_10m_fl4_L49
## 3                             StandardizedHeight_10m          elev_ws5_10m
## 4                            geom_dtm_10m_hyd_fl5_L5   Catchment_slope_50m
## 5                                    geom_10m_fl1_L4      slope_ws29_hr_hr
## 6    geom_hr_L50m_fl10_rpliprichness_UE_hr_5cells_hr FlowLineCurvature_10m
## [1] "OBB error with all predictors of  allpreds is  0.518518518518518"
##                                                MeanDecreaseGini
## Normalized_Height_50m                                 0.8388551
## Modified_Catchment_Area_50m                           0.6094262
## Standardized_Height_50m                               0.6575744
## Vertical_Distance_to_Channel_Network_50m              0.3790405
## Catchment_Area2_50m                                   0.3992647
## SGUT_wTGnew                                           0.4280335
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1        0.2957347
## vectorstrength_hr_ws33_hr                             0.1977812
## Topographic_Wetness_Index_50m                         0.2297360
## SGUT_wTGnew.1                                         0.4821263
##                                                MeanDecreaseAccuracy
## Normalized_Height_50m                                   0.005078172
## Modified_Catchment_Area_50m                             0.002624196
## Standardized_Height_50m                                 0.002144330
## Vertical_Distance_to_Channel_Network_50m                0.001967351
## Catchment_Area2_50m                                     0.001936848
## SGUT_wTGnew                                             0.001698795
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1          0.001624830
## vectorstrength_hr_ws33_hr                               0.001556860
## Topographic_Wetness_Index_50m                           0.001488524
## SGUT_wTGnew.1                                           0.001394048
##         
## altpreds  1  2  3  4  5
##        1 39  0  0  0  0
##        2  0 24  0  0  0
##        3  0  0 20  0  0
##        4  0  0  0 16  0
##        5  0  0  0  0  9
## [1] "classification error rate with altdata:  0"
## [1] "OBB error with all predictors of  allpreds is  0.546296296296296"
## [1] "OOB-error:  0.537037037037037  for predictors geom_dtm_10m_hyd_fl5_L110 AND Catchment_area_50m"
##      
## preds  1  2  3  4  5
##     1 37  6  4  0  1
##     2  0 15  1  2  0
##     3  2  3 14  1  3
##     4  0  0  1 13  1
##     5  0  0  0  0  4
## [1] "Kappa overall =  0.686920222634508"
## [1] "Tau overall =  0.710648148148148"
## [1] "mean quality =  0.588996168582375"
## [1] "The quality is  0.588996168582375"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.702038641647599"
## [1] "classification error rate with altdata:  0.231481481481481"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  allpreds is  0.583333333333333"
## [1] "OOB-error:  0.416666666666667  for predictors geom_dtm_10m_hyd_fl5_L110 AND Catchment_area_50m AND geom_10m_fl4_L40 AND planc_ws5_hr_hr AND geom_dtm_10m_hyd_fl5_L42 AND terraintexture_hr_ws53_tp25 AND geom_dtm_10m_hyd_fl5_L140 AND planc_DTM_50m_avg_ws9_50m AND SGU AND crosc_ws5_hr_hr AND geom_10m_fl10_L3 AND geom_dtm_10m_hyd_fl5_L4 AND geom_hr_L3_fl10_rplipshannon_UE_hr_10cells_hr AND geom_dtm_10m_hyd_fl5_L110 AND maxic_ws3_hr_hr AND terraintexture_hr_ws41_tp5 AND planc_ws29_hr_hr AND LS_Factor_50m AND geom_10m_fl1_L12 AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr AND StandardizedHeight_10m AND geom_dtm_10m_hyd_fl5_L5 AND geom_10m_fl1_L4 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_5cells_hr AND geom_10m_fl1_L8 AND geom_10m_fl4_L49 AND elev_ws5_10m AND Catchment_slope_50m AND slope_ws29_hr_hr AND FlowLineCurvature_10m"
##      
## preds  1  2  3  4  5
##     1 39  0  0  0  0
##     2  0 24  0  0  0
##     3  0  0 20  0  0
##     4  0  0  0 16  0
##     5  0  0  0  0  9
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
i=1
for(i in 1:length(dependentlist)){
  dep=dependentlist[i]
  print(dependentlist_eng[i])
  load(paste("../data/modeldata/RForigmodeldatawithgeoandgeom_",dep,".RData",sep=""))  
  preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/RF_localterrain/RFwithgeoandgeom_fw_5fold_6p_",dep,"_localterrain",sep=""),kk=1:5,endround = 6,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
  print(preds)
  importance_ranfor_pset(modeldata=origmodeldata,dependent=dep,pset=1,altdata=origmodeldata,withalt=F)
predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
  predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=1,kappasum = F,tausum = F,altdata=origmodeldata)
  predictors = vector()
  for (n in 1:ncol(preds)){
    for (n2 in 1:nrow(preds)){
      predictors <- c(predictors,as.character(preds[n2,n]))
    }
  } 
    print("###################### WITH  PREDICTORS from the FW SELECTION ###################")
  uniquepredictors <- unique(predictors)
  predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=1,kappasum = F,tausum = F,altdata=origmodeldata)
  
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"
## [1] "Prediction error at end is:  0.703463203463203"
## [2] "Prediction error at end is:  0.731168831168831"
## [3] "Prediction error at end is:  0.704329004329004"
## [4] "Prediction error at end is:  0.686580086580087"
## [5] "Prediction error at end is:  0.667532467532468"
## [6] "Prediction error at end is:  0.685714285714286"
##                          k 1                        k 2
## 1                  Slope_50m           minic_ws23_hr_hr
## 2 CrossSectionalCurvature_hr           profc_ws13_hr_hr
## 3  longc_DTM_50m_avg_ws3_50m slope_DTM_50m_avg_ws11_50m
## 4       MaximalCurvature_10m              Convexity_50m
## 5  longc_DTM_50m_avg_ws5_50m           minic_ws19_hr_hr
## 6  profc_DTM_50m_avg_ws3_50m           profc_ws11_hr_hr
##                          k 3                       k 4
## 1            maxic_ws9_hr_hr             crosc_ws3_10m
## 2              crosc_ws7_10m       Total_Curvature_50m
## 3 slope_DTM_50m_avg_ws11_50m         Plan_Curvature_hr
## 4              planc_ws5_10m planc_DTM_50m_avg_ws3_50m
## 5            maxic_ws5_hr_hr           planc_ws3_hr_hr
## 6           minic_ws11_hr_hr                 Slope_50m
##                    k 5
## 1     crosc_ws19_hr_hr
## 2    Plan_Curvature_hr
## 3      longc_ws5_hr_hr
## 4     crosc_ws23_hr_hr
## 5      maxic_ws3_hr_hr
## 6 Maximal_Curvature_hr
## [1] "OBB error with all predictors of  localterrain is  0.601851851851852"
##                           MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws3_10m                    1.0265466          0.003632591
## slope_ws5_10m                    0.9122791          0.003233536
## slope_ws13_hr_hr                 0.8457388          0.003166309
## Slope_50m                        0.9580178          0.002916099
## slope_ws15_10m                   0.7347631          0.002781884
## slope_DTM_50m_avg_ws9_50m        0.8000160          0.002553896
## slope_ws11_hr_hr                 0.8156472          0.002542512
## crosc_ws19_hr_hr                 0.5710928          0.002400407
## Slope_10m                        0.8598335          0.002237079
## maxic_ws11_hr_hr                 0.5930079          0.002120002
## [1] "OBB error with all predictors of  localterrain is  0.62037037037037"
## [1] "OOB-error:  0.564814814814815  for predictors Slope_50m AND CrossSectionalCurvature_hr"
##      
## preds  1  2  3  4  5
##     1  4  0  0  0  0
##     2  0 13  0  0  0
##     3  0  0 21  0  0
##     4  0  0  0 38  0
##     5  0  0  0  0 32
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.611111111111111"
## [1] "OOB-error:  0.574074074074074  for predictors Slope_50m AND CrossSectionalCurvature_hr AND longc_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND longc_DTM_50m_avg_ws5_50m AND profc_DTM_50m_avg_ws3_50m AND minic_ws23_hr_hr AND profc_ws13_hr_hr AND slope_DTM_50m_avg_ws11_50m AND Convexity_50m AND minic_ws19_hr_hr AND profc_ws11_hr_hr AND maxic_ws9_hr_hr AND crosc_ws7_10m AND slope_DTM_50m_avg_ws11_50m AND planc_ws5_10m AND maxic_ws5_hr_hr AND minic_ws11_hr_hr AND crosc_ws3_10m AND Total_Curvature_50m AND Plan_Curvature_hr AND planc_DTM_50m_avg_ws3_50m AND planc_ws3_hr_hr AND Slope_50m AND crosc_ws19_hr_hr AND Plan_Curvature_hr AND longc_ws5_hr_hr AND crosc_ws23_hr_hr AND maxic_ws3_hr_hr AND Maximal_Curvature_hr"
##      
## preds  1  2  3  4  5
##     1  4  0  0  0  0
##     2  0 13  0  0  0
##     3  0  0 21  0  0
##     4  0  0  0 38  0
##     5  0  0  0  0 32
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"
## [1] "Prediction error at end is:  0.769264069264069"
## [2] "Prediction error at end is:  0.546320346320346"
## [3] "Prediction error at end is:  0.564502164502165"
## [4] "Prediction error at end is:  0.584415584415584"
## [5] "Prediction error at end is:  0.584415584415584"
## [6] "Prediction error at end is:  0.583982683982684"
##                         k 1                       k 2
## 1                 Slope_50m         Plan_Curvature_hr
## 2      GeneralCurvature_10m          crosc_ws13_hr_hr
## 3    Flow_Line_Curvature_hr          maxic_ws13_hr_hr
## 4 LongitudinalCurvature_10m slope_DTM_50m_avg_ws3_50m
## 5      Minimal_Curvature_hr          slope_ws29_hr_hr
## 6             longc_ws5_10m     Profile_Curvature_50m
##                         k 3                       k 4
## 1       Total_Curvature_50m           maxic_ws9_hr_hr
## 2             Convexity_50m slope_DTM_50m_avg_ws5_50m
## 3            maxic_ws11_10m           maxic_ws5_hr_hr
## 4 slope_DTM_50m_avg_ws3_50m          maxic_ws19_hr_hr
## 5            profc_ws15_10m          maxic_ws13_hr_hr
## 6            profc_ws11_10m             planc_ws7_10m
##                       k 5
## 1        crosc_ws13_hr_hr
## 2        maxic_ws29_hr_hr
## 3         maxic_ws9_hr_hr
## 4         planc_ws9_hr_hr
## 5 Tangential_Curvature_hr
## 6        slope_ws29_hr_hr
## [1] "OBB error with all predictors of  localterrain is  0.509259259259259"
##                           MeanDecreaseGini MeanDecreaseAccuracy
## Slope_50m                        1.1373329          0.006772203
## maxic_ws9_hr_hr                  0.9402189          0.005657172
## Slope_10m                        0.9382462          0.005540153
## slope_ws11_10m                   0.6850574          0.004272931
## slope_DTM_50m_avg_ws9_50m        1.1119233          0.004260792
## slope_ws15_10m                   0.7460380          0.004107919
## maxic_ws11_hr_hr                 0.7214319          0.003908498
## MaximalCurvature_10m             0.9016103          0.003164278
## slope_ws5_10m                    0.7458977          0.003064262
## GeneralCurvature_10m             0.6611461          0.002844019
## [1] "OBB error with all predictors of  localterrain is  0.481481481481481"
## [1] "OOB-error:  0.481481481481481  for predictors Slope_50m AND GeneralCurvature_10m"
##      
## preds  2  3  4  5
##     2 33  0  0  0
##     3  0 23  0  0
##     4  0  0 35  0
##     5  0  0  0 17
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.472222222222222"
## [1] "OOB-error:  0.425925925925926  for predictors Slope_50m AND GeneralCurvature_10m AND Flow_Line_Curvature_hr AND LongitudinalCurvature_10m AND Minimal_Curvature_hr AND longc_ws5_10m AND Plan_Curvature_hr AND crosc_ws13_hr_hr AND maxic_ws13_hr_hr AND slope_DTM_50m_avg_ws3_50m AND slope_ws29_hr_hr AND Profile_Curvature_50m AND Total_Curvature_50m AND Convexity_50m AND maxic_ws11_10m AND slope_DTM_50m_avg_ws3_50m AND profc_ws15_10m AND profc_ws11_10m AND maxic_ws9_hr_hr AND slope_DTM_50m_avg_ws5_50m AND maxic_ws5_hr_hr AND maxic_ws19_hr_hr AND maxic_ws13_hr_hr AND planc_ws7_10m AND crosc_ws13_hr_hr AND maxic_ws29_hr_hr AND maxic_ws9_hr_hr AND planc_ws9_hr_hr AND Tangential_Curvature_hr AND slope_ws29_hr_hr"
##      
## preds  2  3  4  5
##     2 33  0  0  0
##     3  0 23  0  0
##     4  0  0 35  0
##     5  0  0  0 17
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1a.3 - Habitat for soil organisms"

## [1] "Prediction error at end is:  0.675324675324675"
## [2] "Prediction error at end is:  0.694372294372294"
## [3] "Prediction error at end is:  0.648917748917749"
## [4] "Prediction error at end is:  0.602597402597403"
## [5] "Prediction error at end is:  0.528571428571429"
## [6] "Prediction error at end is:  0.567099567099567"
##                         k 1                        k 2
## 1           minic_ws9_hr_hr            minic_ws9_hr_hr
## 2          slope_ws19_hr_hr  profc_DTM_50m_avg_ws3_50m
## 3 minic_DTM_50m_avg_ws7_50m           crosc_ws13_hr_hr
## 4 slope_DTM_50m_avg_ws9_50m          dtm_hr_CONVEX_r30
## 5           minic_ws5_hr_hr           planc_ws15_hr_hr
## 6           planc_ws9_hr_hr crosc_DTM_50m_avg_ws11_50m
##                           k 3                       k 4
## 1   crosc_DTM_50m_avg_ws7_50m     FlowLineCurvature_10m
## 2            crosc_ws15_hr_hr       Total_Curvature_50m
## 3   slope_DTM_50m_avg_ws3_50m            profc_ws15_10m
## 4               Convexity_50m      Convergence_Index_hr
## 5   planc_DTM_50m_avg_ws5_50m slope_DTM_50m_avg_ws7_50m
## 6 CrossSectionalCurvature_50m           maxic_ws3_hr_hr
##                          k 5
## 1            minic_ws9_hr_hr
## 2             slope_ws15_10m
## 3          dtm_hr_CONVEX_r30
## 4 profc_DTM_50m_avg_ws11_50m
## 5           minic_ws23_hr_hr
## 6             crosc_ws15_10m
## [1] "OBB error with all predictors of  localterrain is  0.453703703703704"
##                           MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m               1.6823580          0.007306444
## slope_DTM_50m_avg_ws3_50m        1.2725379          0.005536008
## minic_ws5_hr_hr                  1.1382578          0.004968556
## maxic_ws15_hr_hr                 0.8943811          0.004615553
## maxic_ws23_hr_hr                 0.7709114          0.004033272
## Slope_50m                        0.8929971          0.003954678
## maxic_DTM_50m_avg_ws9_50m        0.8513740          0.003878230
## slope_ws15_10m                   0.9317333          0.003659854
## slope_DTM_50m_avg_ws5_50m        0.9155822          0.003635030
## maxic_DTM_50m_avg_ws5_50m        0.6806564          0.003433799
## [1] "OBB error with all predictors of  localterrain is  0.444444444444444"
## [1] "OOB-error:  0.564814814814815  for predictors minic_ws9_hr_hr AND slope_ws19_hr_hr"
##      
## preds  1  2  3  4
##     1 24  0  0  0
##     2  0 34  0  0
##     3  0  0 26  0
##     4  0  0  0 24
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.490740740740741"
## [1] "OOB-error:  0.444444444444444  for predictors minic_ws9_hr_hr AND slope_ws19_hr_hr AND minic_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws9_50m AND minic_ws5_hr_hr AND planc_ws9_hr_hr AND minic_ws9_hr_hr AND profc_DTM_50m_avg_ws3_50m AND crosc_ws13_hr_hr AND dtm_hr_CONVEX_r30 AND planc_ws15_hr_hr AND crosc_DTM_50m_avg_ws11_50m AND crosc_DTM_50m_avg_ws7_50m AND crosc_ws15_hr_hr AND slope_DTM_50m_avg_ws3_50m AND Convexity_50m AND planc_DTM_50m_avg_ws5_50m AND CrossSectionalCurvature_50m AND FlowLineCurvature_10m AND Total_Curvature_50m AND profc_ws15_10m AND Convergence_Index_hr AND slope_DTM_50m_avg_ws7_50m AND maxic_ws3_hr_hr AND minic_ws9_hr_hr AND slope_ws15_10m AND dtm_hr_CONVEX_r30 AND profc_DTM_50m_avg_ws11_50m AND minic_ws23_hr_hr AND crosc_ws15_10m"
##      
## preds  1  2  3  4
##     1 24  0  0  0
##     2  0 34  0  0
##     3  0  0 26  0
##     4  0  0  0 24
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1a.4 - Habitat for crops"
## [1] "Prediction error at end is:  0.267099567099567"
## [2] "Prediction error at end is:  0.183982683982684"
## [3] "Prediction error at end is:  0.174458874458874"
## [4] "Prediction error at end is:  0.212121212121212"
## [5] "Prediction error at end is:  0.212121212121212"
## [6] "Prediction error at end is:  0.211688311688312"
##                k 1                     k 2                       k 3
## 1  slope_ws5_hr_hr           slope_ws5_10m           slope_ws3_hr_hr
## 2 crosc_ws13_hr_hr        maxic_ws13_hr_hr          longc_ws15_hr_hr
## 3         Slope_hr           slope_ws3_10m minic_DTM_50m_avg_ws9_50m
## 4 minic_ws23_hr_hr          minic_ws11_10m             longc_ws7_10m
## 5    minic_ws5_10m Tangential_Curvature_hr           longc_ws3_hr_hr
## 6 planc_ws19_hr_hr         maxic_ws7_hr_hr profc_DTM_50m_avg_ws7_50m
##                         k 4                       k 5
## 1           slope_ws3_hr_hr                  Slope_hr
## 2          minic_ws15_hr_hr          longc_ws15_hr_hr
## 3 minic_DTM_50m_avg_ws9_50m           slope_ws5_hr_hr
## 4             Convexity_10m          profc_ws15_hr_hr
## 5           longc_ws7_hr_hr slope_DTM_50m_avg_ws9_50m
## 6          minic_ws23_hr_hr      ProfileCurvature_10m
## [1] "OBB error with all predictors of  localterrain is  0.157407407407407"
##                  MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws3_hr_hr          3.090598           0.02574590
## Slope_hr                 3.030718           0.02508142
## slope_ws5_hr_hr          2.709475           0.02110797
## slope_ws23_hr_hr         2.400080           0.01759921
## Slope_10m                2.016602           0.01630050
## slope_ws5_10m            1.773646           0.01471760
## slope_ws19_hr_hr         1.977858           0.01424233
## slope_ws13_hr_hr         2.081435           0.01402077
## slope_ws11_hr_hr         2.180392           0.01332207
## slope_ws7_10m            1.696628           0.01207609
## [1] "OBB error with all predictors of  localterrain is  0.185185185185185"
## [1] "OOB-error:  0.12037037037037  for predictors slope_ws5_hr_hr AND crosc_ws13_hr_hr"
##      
## preds  3  4  5
##     3  8  0  0
##     4  0 62  0
##     5  0  0 38
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.175925925925926"
## [1] "OOB-error:  0.12037037037037  for predictors slope_ws5_hr_hr AND crosc_ws13_hr_hr AND Slope_hr AND minic_ws23_hr_hr AND minic_ws5_10m AND planc_ws19_hr_hr AND slope_ws5_10m AND maxic_ws13_hr_hr AND slope_ws3_10m AND minic_ws11_10m AND Tangential_Curvature_hr AND maxic_ws7_hr_hr AND slope_ws3_hr_hr AND longc_ws15_hr_hr AND minic_DTM_50m_avg_ws9_50m AND longc_ws7_10m AND longc_ws3_hr_hr AND profc_DTM_50m_avg_ws7_50m AND slope_ws3_hr_hr AND minic_ws15_hr_hr AND minic_DTM_50m_avg_ws9_50m AND Convexity_10m AND longc_ws7_hr_hr AND minic_ws23_hr_hr AND Slope_hr AND longc_ws15_hr_hr AND slope_ws5_hr_hr AND profc_ws15_hr_hr AND slope_DTM_50m_avg_ws9_50m AND ProfileCurvature_10m"
##      
## preds  3  4  5
##     3  8  0  0
##     4  0 62  0
##     5  0  0 38
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.1 - Average precipitation retention capacity"
## [1] "Prediction error at end is:  0.751082251082251"
## [2] "Prediction error at end is:  0.674458874458874"
## [3] "Prediction error at end is:  0.629437229437229"
## [4] "Prediction error at end is:  0.583982683982684"
## [5] "Prediction error at end is:  0.601298701298701"
## [6] "Prediction error at end is:  0.592207792207792"
##                          k 1                       k 2
## 1           crosc_ws23_hr_hr crosc_DTM_50m_avg_ws9_50m
## 2              planc_ws5_10m             profc_ws7_10m
## 3  maxic_DTM_50m_avg_ws3_50m             crosc_ws5_10m
## 4 slope_DTM_50m_avg_ws11_50m      Convergence_Index_hr
## 5  slope_DTM_50m_avg_ws5_50m planc_DTM_50m_avg_ws9_50m
## 6          dtm_hr_CONVEX_r10            slope_ws15_10m
##                          k 3              k 4                       k 5
## 1           slope_ws29_hr_hr crosc_ws23_hr_hr           slope_ws3_hr_hr
## 2            slope_ws5_hr_hr slope_ws11_hr_hr             maxic_ws3_10m
## 3       Minimal_Curvature_hr    maxic_ws3_10m slope_DTM_50m_avg_ws9_50m
## 4 crosc_DTM_50m_avg_ws11_50m maxic_ws23_hr_hr           maxic_ws5_hr_hr
## 5             minic_ws15_10m planc_ws19_hr_hr      ProfileCurvature_10m
## 6         TotalCurvature_10m        Slope_50m            slope_ws15_10m
## [1] "OBB error with all predictors of  localterrain is  0.648148148148148"
##                            MeanDecreaseGini MeanDecreaseAccuracy
## MaximalCurvature_10m              0.7707003          0.002423896
## profc_ws3_hr_hr                   0.7067442          0.002401000
## Minimal_Curvature_hr              0.5479792          0.002055238
## Slope_hr                          0.6140058          0.002037640
## slope_DTM_50m_avg_ws11_50m        0.7824028          0.001692910
## Flow_Line_Curvature_50m           0.7049653          0.001623029
## longc_ws3_hr_hr                   0.7170861          0.001527864
## ProfileCurvature_10m              0.6099676          0.001442390
## maxic_ws5_10m                     0.5446579          0.001433610
## minic_ws5_hr_hr                   0.8674175          0.001372780
## [1] "OBB error with all predictors of  localterrain is  0.62037037037037"
## [1] "OOB-error:  0.546296296296296  for predictors crosc_ws23_hr_hr AND planc_ws5_10m"
##      
## preds  1  2  3  4  5
##     1 26  0  0  0  0
##     2  0 45  0  0  0
##     3  0  0 12  0  0
##     4  0  0  0 18  0
##     5  0  0  0  0  7
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.62037037037037"
## [1] "OOB-error:  0.574074074074074  for predictors crosc_ws23_hr_hr AND planc_ws5_10m AND maxic_DTM_50m_avg_ws3_50m AND slope_DTM_50m_avg_ws11_50m AND slope_DTM_50m_avg_ws5_50m AND dtm_hr_CONVEX_r10 AND crosc_DTM_50m_avg_ws9_50m AND profc_ws7_10m AND crosc_ws5_10m AND Convergence_Index_hr AND planc_DTM_50m_avg_ws9_50m AND slope_ws15_10m AND slope_ws29_hr_hr AND slope_ws5_hr_hr AND Minimal_Curvature_hr AND crosc_DTM_50m_avg_ws11_50m AND minic_ws15_10m AND TotalCurvature_10m AND crosc_ws23_hr_hr AND slope_ws11_hr_hr AND maxic_ws3_10m AND maxic_ws23_hr_hr AND planc_ws19_hr_hr AND Slope_50m AND slope_ws3_hr_hr AND maxic_ws3_10m AND slope_DTM_50m_avg_ws9_50m AND maxic_ws5_hr_hr AND ProfileCurvature_10m AND slope_ws15_10m"
##      
## preds  1  2  3  4  5
##     1 26  0  0  0  0
##     2  0 45  0  0  0
##     3  0  0 12  0  0
##     4  0  0  0 18  0
##     5  0  0  0  0  7
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.1 - Minimum precipitation retention capacity"
## [1] "Prediction error at end is:  0.814718614718615"
## [2] "Prediction error at end is:  0.693073593073593"
## [3] "Prediction error at end is:  0.75021645021645" 
## [4] "Prediction error at end is:  0.694372294372294"
## [5] "Prediction error at end is:  0.675324675324675"
## [6] "Prediction error at end is:  0.702597402597403"
##                  k 1                     k 2              k 3
## 1      profc_ws7_10m        longc_ws23_hr_hr profc_ws29_hr_hr
## 2    crosc_ws7_hr_hr        crosc_ws29_hr_hr maxic_ws15_hr_hr
## 3   planc_ws29_hr_hr        minic_ws15_hr_hr    minic_ws7_10m
## 4 Total_Curvature_hr           crosc_ws5_10m    profc_ws3_10m
## 5      crosc_ws3_10m Flow_Line_Curvature_50m minic_ws13_hr_hr
## 6   minic_ws15_hr_hr          profc_ws11_10m maxic_ws29_hr_hr
##                          k 4                        k 5
## 1              maxic_ws7_10m           planc_ws29_hr_hr
## 2           planc_ws13_hr_hr              longc_ws5_10m
## 3              crosc_ws5_10m           crosc_ws15_hr_hr
## 4  longc_DTM_50m_avg_ws9_50m             crosc_ws15_10m
## 5 CrossSectionalCurvature_hr planc_DTM_50m_avg_ws11_50m
## 6   Tangential_Curvature_50m              profc_ws5_10m
## [1] "OBB error with all predictors of  localterrain is  0.666666666666667"
##                       MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m           1.6915471          0.004384354
## crosc_ws5_10m                0.6672548          0.003256847
## maxic_ws13_hr_hr             0.6509108          0.002767572
## Total_Curvature_50m          1.2497922          0.002696361
## crosc_ws11_10m               0.6882563          0.002473567
## crosc_ws7_10m                0.4615653          0.002408718
## maxic_ws15_hr_hr             0.7767013          0.001615967
## FlowLineCurvature_10m        0.6296937          0.001611732
## maxic_ws11_hr_hr             0.6001174          0.001608941
## planc_ws11_10m               0.5128324          0.001540164
## [1] "OBB error with all predictors of  localterrain is  0.62962962962963"
## [1] "OOB-error:  0.657407407407407  for predictors profc_ws7_10m AND crosc_ws7_hr_hr"
##      
## preds  1  2  3  4  5
##     1 12  0  0  0  0
##     2  0 26  0  0  0
##     3  0  0 22  0  0
##     4  0  0  0 26  0
##     5  0  0  0  0 22
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.666666666666667"
## [1] "OOB-error:  0.490740740740741  for predictors profc_ws7_10m AND crosc_ws7_hr_hr AND planc_ws29_hr_hr AND Total_Curvature_hr AND crosc_ws3_10m AND minic_ws15_hr_hr AND longc_ws23_hr_hr AND crosc_ws29_hr_hr AND minic_ws15_hr_hr AND crosc_ws5_10m AND Flow_Line_Curvature_50m AND profc_ws11_10m AND profc_ws29_hr_hr AND maxic_ws15_hr_hr AND minic_ws7_10m AND profc_ws3_10m AND minic_ws13_hr_hr AND maxic_ws29_hr_hr AND maxic_ws7_10m AND planc_ws13_hr_hr AND crosc_ws5_10m AND longc_DTM_50m_avg_ws9_50m AND CrossSectionalCurvature_hr AND Tangential_Curvature_50m AND planc_ws29_hr_hr AND longc_ws5_10m AND crosc_ws15_hr_hr AND crosc_ws15_10m AND planc_DTM_50m_avg_ws11_50m AND profc_ws5_10m"
##      
## preds  1  2  3  4  5
##     1 12  0  0  0  0
##     2  0 26  0  0  0
##     3  0  0 22  0  0
##     4  0  0  0 26  0
##     5  0  0  0  0 22
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.2 - Retention capacity for heavy precipitation events"
## [1] "Prediction error at end is:  0.471428571428571"
## [2] "Prediction error at end is:  0.352380952380952"
## [3] "Prediction error at end is:  0.287878787878788"
## [4] "Prediction error at end is:  0.296536796536797"
## [5] "Prediction error at end is:  0.296536796536797"
## [6] "Prediction error at end is:  0.296969696969697"
##                         k 1                       k 2
## 1     Minimal_Curvature_50m          maxic_ws23_hr_hr
## 2 Longitudinal_Curvature_hr          minic_ws29_hr_hr
## 3 longc_DTM_50m_avg_ws3_50m             Convexity_50m
## 4      MaximalCurvature_10m slope_DTM_50m_avg_ws9_50m
## 5          maxic_ws13_hr_hr         Plan_Curvature_hr
## 6            minic_ws15_10m      MaximalCurvature_10m
##                         k 3                    k 4
## 1             crosc_ws7_10m          Convexity_50m
## 2 slope_DTM_50m_avg_ws9_50m         longc_ws15_10m
## 3           crosc_ws5_hr_hr  FlowLineCurvature_10m
## 4           minic_ws5_hr_hr        minic_ws3_hr_hr
## 5 profc_DTM_50m_avg_ws5_50m         profc_ws11_10m
## 6        Plan_Curvature_50m Flow_Line_Curvature_hr
##                          k 5
## 1           maxic_ws23_hr_hr
## 2           minic_ws23_hr_hr
## 3      Minimal_Curvature_50m
## 4 maxic_DTM_50m_avg_ws11_50m
## 5 profc_DTM_50m_avg_ws11_50m
## 6             minic_ws15_10m
## [1] "OBB error with all predictors of  localterrain is  0.277777777777778"
##                           MeanDecreaseGini MeanDecreaseAccuracy
## FlowLineCurvature_10m            0.9386443          0.006205879
## crosc_DTM_50m_avg_ws5_50m        0.8284482          0.005603521
## minic_ws7_10m                    0.7772594          0.005294088
## minic_ws29_hr_hr                 0.8524506          0.004889761
## crosc_ws11_10m                   0.6672203          0.003382759
## maxic_ws7_10m                    0.3935556          0.003287878
## maxic_ws23_hr_hr                 0.3907311          0.003249979
## minic_ws23_hr_hr                 0.6889049          0.003240223
## minic_ws19_hr_hr                 0.5945955          0.003059828
## longc_ws5_10m                    0.5116574          0.002863669
## [1] "OBB error with all predictors of  localterrain is  0.305555555555556"
## [1] "OOB-error:  0.240740740740741  for predictors Minimal_Curvature_50m AND Longitudinal_Curvature_hr"
##      
## preds  1  2  3  4  5
##     1 75  0  0  0  0
##     2  0 11  0  0  0
##     3  0  0  6  0  0
##     4  0  0  0  8  0
##     5  0  0  0  0  8
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.277777777777778"
## [1] "OOB-error:  0.25  for predictors Minimal_Curvature_50m AND Longitudinal_Curvature_hr AND longc_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND maxic_ws13_hr_hr AND minic_ws15_10m AND maxic_ws23_hr_hr AND minic_ws29_hr_hr AND Convexity_50m AND slope_DTM_50m_avg_ws9_50m AND Plan_Curvature_hr AND MaximalCurvature_10m AND crosc_ws7_10m AND slope_DTM_50m_avg_ws9_50m AND crosc_ws5_hr_hr AND minic_ws5_hr_hr AND profc_DTM_50m_avg_ws5_50m AND Plan_Curvature_50m AND Convexity_50m AND longc_ws15_10m AND FlowLineCurvature_10m AND minic_ws3_hr_hr AND profc_ws11_10m AND Flow_Line_Curvature_hr AND maxic_ws23_hr_hr AND minic_ws23_hr_hr AND Minimal_Curvature_50m AND maxic_DTM_50m_avg_ws11_50m AND profc_DTM_50m_avg_ws11_50m AND minic_ws15_10m"
##      
## preds  1  2  3  4  5
##     1 75  0  0  0  0
##     2  0 11  0  0  0
##     3  0  0  6  0  0
##     4  0  0  0  8  0
##     5  0  0  0  0  8
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.3 - groundwater reformation rate"
## [1] "Prediction error at end is:  0.788311688311688"
## [2] "Prediction error at end is:  0.73030303030303" 
## [3] "Prediction error at end is:  0.693506493506494"
## [4] "Prediction error at end is:  0.700865800865801"
## [5] "Prediction error at end is:  0.674025974025974"
## [6] "Prediction error at end is:  0.701731601731602"
##                         k 1                       k 2
## 1          slope_ws15_hr_hr LongitudinalCurvature_10m
## 2          profc_ws23_hr_hr             maxic_ws7_10m
## 3          maxic_ws19_hr_hr           maxic_ws3_hr_hr
## 4 slope_DTM_50m_avg_ws5_50m          maxic_ws29_hr_hr
## 5 longc_DTM_50m_avg_ws3_50m slope_DTM_50m_avg_ws5_50m
## 6          minic_ws29_hr_hr      GeneralCurvature_10m
##                         k 3                       k 4
## 1            minic_ws15_10m          minic_ws11_hr_hr
## 2     FlowLineCurvature_10m maxic_DTM_50m_avg_ws5_50m
## 3          dtm_hr_CONVEX_r5          crosc_ws13_hr_hr
## 4            planc_ws11_10m slope_DTM_50m_avg_ws9_50m
## 5 minic_DTM_50m_avg_ws3_50m             planc_ws3_10m
## 6             crosc_ws7_10m           maxic_ws9_hr_hr
##                         k 5
## 1            minic_ws15_10m
## 2          maxic_ws19_hr_hr
## 3           longc_ws3_hr_hr
## 4          maxic_ws15_hr_hr
## 5 maxic_DTM_50m_avg_ws5_50m
## 6             planc_ws7_10m
## [1] "OBB error with all predictors of  localterrain is  0.611111111111111"
##                      MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws5_hr_hr             0.9091755          0.004505613
## slope_ws23_hr_hr            0.8713569          0.003693752
## slope_ws19_hr_hr            0.6994681          0.003065583
## Profile_Curvature_hr        0.8664606          0.003013326
## slope_ws7_10m               0.7242359          0.002625295
## slope_ws15_hr_hr            0.6548519          0.002503940
## minic_ws3_hr_hr             0.9761192          0.002479712
## slope_ws5_10m               0.8099342          0.002418275
## Slope_10m                   0.7466590          0.002312560
## Minimal_Curvature_hr        0.5855733          0.002036061
## [1] "OBB error with all predictors of  localterrain is  0.611111111111111"
## [1] "OOB-error:  0.62037037037037  for predictors slope_ws15_hr_hr AND profc_ws23_hr_hr"
##      
## preds  1  2  3  4  5
##     1 14  0  0  0  0
##     2  0 39  0  0  0
##     3  0  0 14  0  0
##     4  0  0  0 19  0
##     5  0  0  0  0 22
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.611111111111111"
## [1] "OOB-error:  0.564814814814815  for predictors slope_ws15_hr_hr AND profc_ws23_hr_hr AND maxic_ws19_hr_hr AND slope_DTM_50m_avg_ws5_50m AND longc_DTM_50m_avg_ws3_50m AND minic_ws29_hr_hr AND LongitudinalCurvature_10m AND maxic_ws7_10m AND maxic_ws3_hr_hr AND maxic_ws29_hr_hr AND slope_DTM_50m_avg_ws5_50m AND GeneralCurvature_10m AND minic_ws15_10m AND FlowLineCurvature_10m AND dtm_hr_CONVEX_r5 AND planc_ws11_10m AND minic_DTM_50m_avg_ws3_50m AND crosc_ws7_10m AND minic_ws11_hr_hr AND maxic_DTM_50m_avg_ws5_50m AND crosc_ws13_hr_hr AND slope_DTM_50m_avg_ws9_50m AND planc_ws3_10m AND maxic_ws9_hr_hr AND minic_ws15_10m AND maxic_ws19_hr_hr AND longc_ws3_hr_hr AND maxic_ws15_hr_hr AND maxic_DTM_50m_avg_ws5_50m AND planc_ws7_10m"
##      
## preds  1  2  3  4  5
##     1 14  0  0  0  0
##     2  0 39  0  0  0
##     3  0  0 14  0  0
##     4  0  0  0 19  0
##     5  0  0  0  0 22
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.4 - Potential for providing nutrients for plants"
## [1] "Prediction error at end is:  0.529004329004329"
## [2] "Prediction error at end is:  0.371428571428571"
## [3] "Prediction error at end is:  0.325108225108225"
## [4] "Prediction error at end is:  0.324675324675325"
## [5] "Prediction error at end is:  0.306493506493506"
## [6] "Prediction error at end is:  0.288311688311688"
##                         k 1                        k 2
## 1           minic_ws5_hr_hr            crosc_ws5_hr_hr
## 2 longc_DTM_50m_avg_ws5_50m profc_DTM_50m_avg_ws11_50m
## 3          maxic_ws15_hr_hr planc_DTM_50m_avg_ws11_50m
## 4          profc_ws29_hr_hr           minic_ws19_hr_hr
## 5            maxic_ws15_10m         Plan_Curvature_50m
## 6            maxic_ws11_10m            minic_ws7_hr_hr
##                           k 3                         k 4
## 1            longc_ws29_hr_hr     TangentialCurvature_10m
## 2               minic_ws7_10m              slope_ws11_10m
## 3               longc_ws7_10m   minic_DTM_50m_avg_ws7_50m
## 4             crosc_ws5_hr_hr   slope_DTM_50m_avg_ws7_50m
## 5            profc_ws29_hr_hr CrossSectionalCurvature_50m
## 6 CrossSectionalCurvature_50m   maxic_DTM_50m_avg_ws5_50m
##                          k 5
## 1           longc_ws19_hr_hr
## 2        Total_Curvature_50m
## 3  minic_DTM_50m_avg_ws7_50m
## 4          dtm_hr_CONVEX_r30
## 5 planc_DTM_50m_avg_ws11_50m
## 6           minic_ws15_hr_hr
## [1] "OBB error with all predictors of  localterrain is  0.305555555555556"
##                         MeanDecreaseGini MeanDecreaseAccuracy
## minic_ws11_hr_hr               1.3896970          0.011297624
## minic_ws9_hr_hr                1.3926995          0.007597971
## crosc_ws5_hr_hr                1.1285306          0.007037121
## minic_ws7_hr_hr                1.5603885          0.006705559
## TotalCurvature_10m             1.2511444          0.006160002
## minic_ws5_hr_hr                0.9093794          0.005050313
## planc_ws5_hr_hr                0.6506936          0.003306029
## TangentialCurvature_10m        0.3848811          0.003082114
## MaximalCurvature_10m           0.4665715          0.003028937
## maxic_ws7_10m                  0.7210332          0.003016545
## [1] "OBB error with all predictors of  localterrain is  0.305555555555556"
## [1] "OOB-error:  0.296296296296296  for predictors minic_ws5_hr_hr AND longc_DTM_50m_avg_ws5_50m"
##      
## preds  1  3  5
##     1 65  0  0
##     3  0 11  0
##     5  0  0 32
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.305555555555556"
## [1] "OOB-error:  0.25  for predictors minic_ws5_hr_hr AND longc_DTM_50m_avg_ws5_50m AND maxic_ws15_hr_hr AND profc_ws29_hr_hr AND maxic_ws15_10m AND maxic_ws11_10m AND crosc_ws5_hr_hr AND profc_DTM_50m_avg_ws11_50m AND planc_DTM_50m_avg_ws11_50m AND minic_ws19_hr_hr AND Plan_Curvature_50m AND minic_ws7_hr_hr AND longc_ws29_hr_hr AND minic_ws7_10m AND longc_ws7_10m AND crosc_ws5_hr_hr AND profc_ws29_hr_hr AND CrossSectionalCurvature_50m AND TangentialCurvature_10m AND slope_ws11_10m AND minic_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws7_50m AND CrossSectionalCurvature_50m AND maxic_DTM_50m_avg_ws5_50m AND longc_ws19_hr_hr AND Total_Curvature_50m AND minic_DTM_50m_avg_ws7_50m AND dtm_hr_CONVEX_r30 AND planc_DTM_50m_avg_ws11_50m AND minic_ws15_hr_hr"
##      
## preds  1  3  5
##     1 65  0  0
##     3  0 11  0
##     5  0  0 32
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1c.5 - Potential as a CO2 sink"
## [1] "Prediction error at end is:  0.522943722943723"
## [2] "Prediction error at end is:  0.486580086580087"
## [3] "Prediction error at end is:  0.421212121212121"
## [4] "Prediction error at end is:  0.42987012987013" 
## [5] "Prediction error at end is:  0.401731601731602"
## [6] "Prediction error at end is:  0.354112554112554"
##                  k 1                        k 2                        k 3
## 1   slope_ws19_hr_hr           slope_ws19_hr_hr CrossSectionalCurvature_hr
## 2  PlanCurvature_10m minic_DTM_50m_avg_ws11_50m         TotalCurvature_10m
## 3    crosc_ws5_hr_hr              planc_ws5_10m      FlowLineCurvature_10m
## 4    maxic_ws3_hr_hr           planc_ws13_hr_hr  planc_DTM_50m_avg_ws7_50m
## 5 TotalCurvature_10m   Tangential_Curvature_50m        Total_Curvature_50m
## 6      planc_ws3_10m           crosc_ws13_hr_hr  profc_DTM_50m_avg_ws5_50m
##                 k 4                      k 5
## 1   crosc_ws5_hr_hr         crosc_ws19_hr_hr
## 2  slope_ws13_hr_hr            profc_ws7_10m
## 3  minic_ws15_hr_hr Tangential_Curvature_50m
## 4     longc_ws7_10m           maxic_ws15_10m
## 5  slope_ws15_hr_hr         maxic_ws29_hr_hr
## 6 PlanCurvature_10m  TangentialCurvature_10m
## [1] "OBB error with all predictors of  localterrain is  0.351851851851852"
##                            MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m                1.8085695          0.011117328
## maxic_ws15_hr_hr                  1.6261493          0.009747365
## minic_ws11_hr_hr                  0.9883106          0.006768397
## MaximalCurvature_10m              1.0611950          0.006696107
## minic_ws3_hr_hr                   0.9256518          0.005123103
## minic_ws5_hr_hr                   1.1148332          0.004956734
## CrossSectionalCurvature_hr        0.7767361          0.004492110
## slope_ws5_10m                     0.6749511          0.004464976
## maxic_ws23_hr_hr                  0.6333793          0.004412928
## Total_Curvature_hr                1.0246174          0.004146658
## [1] "OBB error with all predictors of  localterrain is  0.37037037037037"
## [1] "OOB-error:  0.37037037037037  for predictors slope_ws19_hr_hr AND PlanCurvature_10m"
##      
## preds  1  2  3  4  5
##     1 51  0  0  0  0
##     2  0  1  0  0  0
##     3  0  0 18  0  0
##     4  0  0  0 33  0
##     5  0  0  0  0  5
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.324074074074074"
## [1] "OOB-error:  0.314814814814815  for predictors slope_ws19_hr_hr AND PlanCurvature_10m AND crosc_ws5_hr_hr AND maxic_ws3_hr_hr AND TotalCurvature_10m AND planc_ws3_10m AND slope_ws19_hr_hr AND minic_DTM_50m_avg_ws11_50m AND planc_ws5_10m AND planc_ws13_hr_hr AND Tangential_Curvature_50m AND crosc_ws13_hr_hr AND CrossSectionalCurvature_hr AND TotalCurvature_10m AND FlowLineCurvature_10m AND planc_DTM_50m_avg_ws7_50m AND Total_Curvature_50m AND profc_DTM_50m_avg_ws5_50m AND crosc_ws5_hr_hr AND slope_ws13_hr_hr AND minic_ws15_hr_hr AND longc_ws7_10m AND slope_ws15_hr_hr AND PlanCurvature_10m AND crosc_ws19_hr_hr AND profc_ws7_10m AND Tangential_Curvature_50m AND maxic_ws15_10m AND maxic_ws29_hr_hr AND TangentialCurvature_10m"
##      
## preds  1  2  3  4  5
##     1 51  0  0  0  0
##     2  0  1  0  0  0
##     3  0  0 18  0  0
##     4  0  0  0 33  0
##     5  0  0  0  0  5
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.1 - Potential for retention of heavy metals"
## [1] "Prediction error at end is:  0.638095238095238"
## [2] "Prediction error at end is:  0.5004329004329"  
## [3] "Prediction error at end is:  0.435930735930736"
## [4] "Prediction error at end is:  0.427705627705628"
## [5] "Prediction error at end is:  0.408658008658009"
## [6] "Prediction error at end is:  0.399134199134199"
##                         k 1                        k 2
## 1            slope_ws15_10m           longc_ws23_hr_hr
## 2             minic_ws7_10m           slope_ws29_hr_hr
## 3         dtm_hr_CONVEX_r30              maxic_ws7_10m
## 4 longc_DTM_50m_avg_ws5_50m CrossSectionalCurvature_hr
## 5 planc_DTM_50m_avg_ws5_50m            crosc_ws5_hr_hr
## 6          crosc_ws15_hr_hr           crosc_ws13_hr_hr
##                       k 3                        k 4
## 1        minic_ws29_hr_hr        Total_Curvature_50m
## 2 Flow_Line_Curvature_50m CrossSectionalCurvature_hr
## 3         maxic_ws5_hr_hr            minic_ws3_hr_hr
## 4         planc_ws3_hr_hr          dtm_hr_CONVEX_r30
## 5       dtm_hr_CONVEX_r20  crosc_DTM_50m_avg_ws7_50m
## 6   FlowLineCurvature_10m         TotalCurvature_10m
##                         k 5
## 1 slope_DTM_50m_avg_ws3_50m
## 2          maxic_ws29_hr_hr
## 3          slope_ws23_hr_hr
## 4        Total_Curvature_hr
## 5          crosc_ws13_hr_hr
## 6      Minimal_Curvature_hr
## [1] "OBB error with all predictors of  localterrain is  0.361111111111111"
##                    MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m        1.4197148          0.007170865
## slope_ws3_10m             0.9638191          0.005567992
## minic_ws5_hr_hr           1.3393282          0.005490887
## maxic_ws15_hr_hr          0.9104475          0.005090625
## Slope_10m                 1.1650392          0.005079061
## slope_ws5_10m             0.8172995          0.004812107
## Slope_50m                 0.6501166          0.004447252
## slope_ws13_hr_hr          0.8421275          0.004091610
## slope_ws23_hr_hr          0.8940395          0.004091067
## crosc_ws13_hr_hr          0.6049590          0.003628726
## [1] "OBB error with all predictors of  localterrain is  0.37037037037037"
## [1] "OOB-error:  0.37962962962963  for predictors slope_ws15_10m AND minic_ws7_10m"
##      
## preds  1  2  3  4  5
##     1 42  0  0  0  0
##     2  0  7  0  0  0
##     3  0  0 10  0  0
##     4  0  0  0  7  0
##     5  0  0  0  0 42
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.37037037037037"
## [1] "OOB-error:  0.333333333333333  for predictors slope_ws15_10m AND minic_ws7_10m AND dtm_hr_CONVEX_r30 AND longc_DTM_50m_avg_ws5_50m AND planc_DTM_50m_avg_ws5_50m AND crosc_ws15_hr_hr AND longc_ws23_hr_hr AND slope_ws29_hr_hr AND maxic_ws7_10m AND CrossSectionalCurvature_hr AND crosc_ws5_hr_hr AND crosc_ws13_hr_hr AND minic_ws29_hr_hr AND Flow_Line_Curvature_50m AND maxic_ws5_hr_hr AND planc_ws3_hr_hr AND dtm_hr_CONVEX_r20 AND FlowLineCurvature_10m AND Total_Curvature_50m AND CrossSectionalCurvature_hr AND minic_ws3_hr_hr AND dtm_hr_CONVEX_r30 AND crosc_DTM_50m_avg_ws7_50m AND TotalCurvature_10m AND slope_DTM_50m_avg_ws3_50m AND maxic_ws29_hr_hr AND slope_ws23_hr_hr AND Total_Curvature_hr AND crosc_ws13_hr_hr AND Minimal_Curvature_hr"
##      
## preds  1  2  3  4  5
##     1 42  0  0  0  0
##     2  0  7  0  0  0
##     3  0  0 10  0  0
##     4  0  0  0  7  0
##     5  0  0  0  0 42
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.2 - Potential for transforming organic contaminants"
## [1] "Prediction error at end is:  0.657142857142857"
## [2] "Prediction error at end is:  0.647186147186147"
## [3] "Prediction error at end is:  0.593073593073593"
## [4] "Prediction error at end is:  0.62034632034632" 
## [5] "Prediction error at end is:  0.603463203463203"
## [6] "Prediction error at end is:  0.601298701298701"
##                         k 1                  k 2                       k 3
## 1          crosc_ws19_hr_hr MaximalCurvature_10m          crosc_ws15_hr_hr
## 2      Minimal_Curvature_hr      maxic_ws3_hr_hr          crosc_ws19_hr_hr
## 3          maxic_ws29_hr_hr MinimalCurvature_10m          maxic_ws15_hr_hr
## 4 minic_DTM_50m_avg_ws9_50m       longc_ws15_10m minic_DTM_50m_avg_ws5_50m
## 5           crosc_ws3_hr_hr     dtm_hr_CONVEX_r5 longc_DTM_50m_avg_ws9_50m
## 6 slope_DTM_50m_avg_ws3_50m     planc_ws15_hr_hr         Plan_Curvature_hr
##                k 4                         k 5
## 1    maxic_ws7_10m            crosc_ws13_hr_hr
## 2 maxic_ws23_hr_hr   crosc_DTM_50m_avg_ws3_50m
## 3   longc_ws11_10m   planc_DTM_50m_avg_ws5_50m
## 4  profc_ws7_hr_hr               crosc_ws3_10m
## 5 minic_ws13_hr_hr CrossSectionalCurvature_50m
## 6    longc_ws7_10m             minic_ws5_hr_hr
## [1] "OBB error with all predictors of  localterrain is  0.638888888888889"
##                      MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m          1.3895773          0.005646626
## MaximalCurvature_10m        1.3814653          0.005120462
## minic_ws5_hr_hr             0.8674868          0.003467020
## slope_ws19_hr_hr            0.7628082          0.003461461
## slope_ws29_hr_hr            0.6657954          0.003205129
## slope_ws5_10m               0.7759054          0.002965819
## minic_ws7_hr_hr             0.8871311          0.002909758
## slope_ws7_10m               0.8566747          0.002885358
## maxic_ws9_hr_hr             0.7843232          0.002630408
## maxic_ws15_hr_hr            0.9504995          0.002571349
## [1] "OBB error with all predictors of  localterrain is  0.601851851851852"
## [1] "OOB-error:  0.555555555555556  for predictors crosc_ws19_hr_hr AND Minimal_Curvature_hr"
##      
## preds  1  2  3  4  5
##     1 10  0  0  0  0
##     2  0 22  0  0  0
##     3  0  0 28  0  0
##     4  0  0  0  9  0
##     5  0  0  0  0 39
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.601851851851852"
## [1] "OOB-error:  0.564814814814815  for predictors crosc_ws19_hr_hr AND Minimal_Curvature_hr AND maxic_ws29_hr_hr AND minic_DTM_50m_avg_ws9_50m AND crosc_ws3_hr_hr AND slope_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND maxic_ws3_hr_hr AND MinimalCurvature_10m AND longc_ws15_10m AND dtm_hr_CONVEX_r5 AND planc_ws15_hr_hr AND crosc_ws15_hr_hr AND crosc_ws19_hr_hr AND maxic_ws15_hr_hr AND minic_DTM_50m_avg_ws5_50m AND longc_DTM_50m_avg_ws9_50m AND Plan_Curvature_hr AND maxic_ws7_10m AND maxic_ws23_hr_hr AND longc_ws11_10m AND profc_ws7_hr_hr AND minic_ws13_hr_hr AND longc_ws7_10m AND crosc_ws13_hr_hr AND crosc_DTM_50m_avg_ws3_50m AND planc_DTM_50m_avg_ws5_50m AND crosc_ws3_10m AND CrossSectionalCurvature_50m AND minic_ws5_hr_hr"
##      
## preds  1  2  3  4  5
##     1 10  0  0  0  0
##     2  0 22  0  0  0
##     3  0  0 28  0  0
##     4  0  0  0  9  0
##     5  0  0  0  0 39
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"
## [1] "Prediction error at end is:  0.0554112554112554"
## [2] "Prediction error at end is:  0.0463203463203463"
## [3] "Prediction error at end is:  0.0372294372294372"
## [4] "Prediction error at end is:  0.0372294372294372"
## [5] "Prediction error at end is:  0.0372294372294372"
## [6] "Prediction error at end is:  0.0372294372294372"
##                     k 1                       k 2
## 1         Convexity_10m         PlanCurvature_10m
## 2         planc_ws5_10m           planc_ws9_hr_hr
## 3         Convexity_50m          slope_ws13_hr_hr
## 4 Maximal_Curvature_50m maxic_DTM_50m_avg_ws9_50m
## 5 Minimal_Curvature_50m slope_DTM_50m_avg_ws3_50m
## 6    Plan_Curvature_50m            slope_ws15_10m
##                           k 3                   k 4
## 1 CrossSectionalCurvature_50m        longc_ws11_10m
## 2           PlanCurvature_10m         longc_ws7_10m
## 3            slope_ws13_hr_hr         profc_ws5_10m
## 4     Flow_Line_Curvature_50m Minimal_Curvature_50m
## 5             maxic_ws7_hr_hr       crosc_ws7_hr_hr
## 6            planc_ws11_hr_hr      maxic_ws23_hr_hr
##                           k 5
## 1               Convexity_10m
## 2             planc_ws9_hr_hr
## 3           PlanCurvature_10m
## 4            minic_ws29_hr_hr
## 5 CrossSectionalCurvature_50m
## 6            minic_ws23_hr_hr
## [1] "OBB error with all predictors of  localterrain is  0.037037037037037"
##                            MeanDecreaseGini MeanDecreaseAccuracy
## profc_DTM_50m_avg_ws5_50m        0.05348294         0.0015212168
## planc_DTM_50m_avg_ws5_50m        0.14353614         0.0014809356
## longc_ws15_10m                   0.05926461         0.0013900200
## planc_DTM_50m_avg_ws11_50m       0.21200915         0.0012811719
## profc_DTM_50m_avg_ws9_50m        0.06087598         0.0011735909
## profc_ws11_10m                   0.05079702         0.0011057148
## profc_ws15_10m                   0.03807431         0.0010371631
## minic_ws29_hr_hr                 0.07496966         0.0009558680
## profc_DTM_50m_avg_ws7_50m        0.10393737         0.0008745406
## longc_DTM_50m_avg_ws7_50m        0.05407145         0.0008603118
## [1] "OBB error with all predictors of  localterrain is  0.037037037037037"
## [1] "OOB-error:  0.037037037037037  for predictors Convexity_10m AND planc_ws5_10m"
##      
## preds   4   5
##     4   4   0
##     5   0 104
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  0.870192307692308"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.037037037037037"
## [1] "OOB-error:  0.0185185185185185  for predictors Convexity_10m AND planc_ws5_10m AND Convexity_50m AND Maximal_Curvature_50m AND Minimal_Curvature_50m AND Plan_Curvature_50m AND PlanCurvature_10m AND planc_ws9_hr_hr AND slope_ws13_hr_hr AND maxic_DTM_50m_avg_ws9_50m AND slope_DTM_50m_avg_ws3_50m AND slope_ws15_10m AND CrossSectionalCurvature_50m AND PlanCurvature_10m AND slope_ws13_hr_hr AND Flow_Line_Curvature_50m AND maxic_ws7_hr_hr AND planc_ws11_hr_hr AND longc_ws11_10m AND longc_ws7_10m AND profc_ws5_10m AND Minimal_Curvature_50m AND crosc_ws7_hr_hr AND maxic_ws23_hr_hr AND Convexity_10m AND planc_ws9_hr_hr AND PlanCurvature_10m AND minic_ws29_hr_hr AND CrossSectionalCurvature_50m AND minic_ws23_hr_hr"
##      
## preds   4   5
##     4   4   0
##     5   0 104
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  0.870192307692308"
## [1] "classification error rate with altdata:  0"
## [1] "1d.4 - Potential for retention of water-soluble contaminants"
## [1] "Prediction error at end is:  0.795238095238095"
## [2] "Prediction error at end is:  0.600865800865801"
## [3] "Prediction error at end is:  0.572727272727273"
## [4] "Prediction error at end is:  0.554545454545455"
## [5] "Prediction error at end is:  0.526406926406926"
## [6] "Prediction error at end is:  0.535930735930736"
##                         k 1                        k 2
## 1      Minimal_Curvature_hr           crosc_ws19_hr_hr
## 2             crosc_ws3_10m CrossSectionalCurvature_hr
## 3           crosc_ws7_hr_hr           crosc_ws23_hr_hr
## 4                 Slope_50m  profc_DTM_50m_avg_ws9_50m
## 5 planc_DTM_50m_avg_ws3_50m  crosc_DTM_50m_avg_ws3_50m
## 6     Profile_Curvature_50m            profc_ws5_hr_hr
##                     k 3                       k 4
## 1  GeneralCurvature_10m          maxic_ws13_hr_hr
## 2 FlowLineCurvature_10m          slope_ws29_hr_hr
## 3       crosc_ws3_hr_hr   Tangential_Curvature_hr
## 4      crosc_ws15_hr_hr LongitudinalCurvature_10m
## 5         longc_ws3_10m          crosc_ws11_hr_hr
## 6  MinimalCurvature_10m longc_DTM_50m_avg_ws3_50m
##                          k 5
## 1           crosc_ws13_hr_hr
## 2       MinimalCurvature_10m
## 3  crosc_DTM_50m_avg_ws9_50m
## 4           crosc_ws15_hr_hr
## 5 crosc_DTM_50m_avg_ws11_50m
## 6            profc_ws7_hr_hr
## [1] "OBB error with all predictors of  localterrain is  0.5"
##                           MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws11_hr_hr                 0.8151606          0.004513121
## slope_ws3_10m                    0.8861537          0.003939477
## TotalCurvature_10m               0.9456823          0.003589144
## slope_ws19_hr_hr                 0.6407110          0.003355717
## slope_ws23_hr_hr                 0.6931134          0.003346353
## slope_DTM_50m_avg_ws5_50m        0.7497832          0.003238108
## Slope_10m                        0.7486609          0.003200858
## crosc_ws7_10m                    0.7558984          0.003080825
## crosc_ws3_10m                    0.5450793          0.002890011
## Minimal_Curvature_hr             0.8398992          0.002520341
## [1] "OBB error with all predictors of  localterrain is  0.5"
## [1] "OOB-error:  0.537037037037037  for predictors Minimal_Curvature_hr AND crosc_ws3_10m"
##      
## preds  1  2  3  4  5
##     1 47  0  0  0  0
##     2  0 14  0  0  0
##     3  0  0 17  0  0
##     4  0  0  0 22  0
##     5  0  0  0  0  8
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.490740740740741"
## [1] "OOB-error:  0.481481481481481  for predictors Minimal_Curvature_hr AND crosc_ws3_10m AND crosc_ws7_hr_hr AND Slope_50m AND planc_DTM_50m_avg_ws3_50m AND Profile_Curvature_50m AND crosc_ws19_hr_hr AND CrossSectionalCurvature_hr AND crosc_ws23_hr_hr AND profc_DTM_50m_avg_ws9_50m AND crosc_DTM_50m_avg_ws3_50m AND profc_ws5_hr_hr AND GeneralCurvature_10m AND FlowLineCurvature_10m AND crosc_ws3_hr_hr AND crosc_ws15_hr_hr AND longc_ws3_10m AND MinimalCurvature_10m AND maxic_ws13_hr_hr AND slope_ws29_hr_hr AND Tangential_Curvature_hr AND LongitudinalCurvature_10m AND crosc_ws11_hr_hr AND longc_DTM_50m_avg_ws3_50m AND crosc_ws13_hr_hr AND MinimalCurvature_10m AND crosc_DTM_50m_avg_ws9_50m AND crosc_ws15_hr_hr AND crosc_DTM_50m_avg_ws11_50m AND profc_ws7_hr_hr"
##      
## preds  1  2  3  4  5
##     1 47  0  0  0  0
##     2  0 14  0  0  0
##     3  0  0 17  0  0
##     4  0  0  0 22  0
##     5  0  0  0  0  8
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "1d.5 - Potential as buffer for acidic contaminants"
## [1] "Prediction error at end is:  0.759307359307359"
## [2] "Prediction error at end is:  0.684848484848485"
## [3] "Prediction error at end is:  0.674458874458874"
## [4] "Prediction error at end is:  0.674891774891775"
## [5] "Prediction error at end is:  0.637662337662338"
## [6] "Prediction error at end is:  0.647186147186147"
##                          k 1              k 2                       k 3
## 1           crosc_ws29_hr_hr crosc_ws23_hr_hr           planc_ws7_hr_hr
## 2      FlowLineCurvature_10m   planc_ws11_10m          maxic_ws23_hr_hr
## 3 planc_DTM_50m_avg_ws11_50m   longc_ws15_10m planc_DTM_50m_avg_ws3_50m
## 4            crosc_ws3_hr_hr    planc_ws3_10m          crosc_ws19_hr_hr
## 5           crosc_ws13_hr_hr   profc_ws15_10m           crosc_ws5_hr_hr
## 6             planc_ws11_10m   crosc_ws15_10m          planc_ws23_hr_hr
##                         k 4              k 5
## 1 minic_DTM_50m_avg_ws9_50m    slope_ws3_10m
## 2           crosc_ws5_hr_hr  planc_ws5_hr_hr
## 3 slope_DTM_50m_avg_ws9_50m crosc_ws19_hr_hr
## 4                 Slope_50m  crosc_ws5_hr_hr
## 5           longc_ws5_hr_hr   maxic_ws11_10m
## 6          longc_ws15_hr_hr slope_ws13_hr_hr
## [1] "OBB error with all predictors of  localterrain is  0.611111111111111"
##                            MeanDecreaseGini MeanDecreaseAccuracy
## maxic_DTM_50m_avg_ws11_50m        1.3473532          0.004675391
## longc_DTM_50m_avg_ws11_50m        0.8828896          0.003151859
## Total_Curvature_50m               0.9422908          0.003136906
## minic_ws11_hr_hr                  0.8122509          0.002817321
## crosc_ws5_hr_hr                   0.9757649          0.002705295
## Slope_10m                         0.8251192          0.002654011
## slope_ws3_10m                     0.6112843          0.002107568
## planc_ws23_hr_hr                  0.6790395          0.002102952
## slope_ws15_hr_hr                  0.6160903          0.002085516
## minic_ws5_10m                     0.5756878          0.002005237
## [1] "OBB error with all predictors of  localterrain is  0.611111111111111"
## [1] "OOB-error:  0.638888888888889  for predictors crosc_ws29_hr_hr AND FlowLineCurvature_10m"
##      
## preds  1  2  3  4  5
##     1 39  0  0  0  0
##     2  0 24  0  0  0
##     3  0  0 20  0  0
##     4  0  0  0 16  0
##     5  0  0  0  0  9
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"
## [1] "###################### WITH  PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of  localterrain is  0.592592592592593"
## [1] "OOB-error:  0.555555555555556  for predictors crosc_ws29_hr_hr AND FlowLineCurvature_10m AND planc_DTM_50m_avg_ws11_50m AND crosc_ws3_hr_hr AND crosc_ws13_hr_hr AND planc_ws11_10m AND crosc_ws23_hr_hr AND planc_ws11_10m AND longc_ws15_10m AND planc_ws3_10m AND profc_ws15_10m AND crosc_ws15_10m AND planc_ws7_hr_hr AND maxic_ws23_hr_hr AND planc_DTM_50m_avg_ws3_50m AND crosc_ws19_hr_hr AND crosc_ws5_hr_hr AND planc_ws23_hr_hr AND minic_DTM_50m_avg_ws9_50m AND crosc_ws5_hr_hr AND slope_DTM_50m_avg_ws9_50m AND Slope_50m AND longc_ws5_hr_hr AND longc_ws15_hr_hr AND slope_ws3_10m AND planc_ws5_hr_hr AND crosc_ws19_hr_hr AND crosc_ws5_hr_hr AND maxic_ws11_10m AND slope_ws13_hr_hr"
##      
## preds  1  2  3  4  5
##     1 39  0  0  0  0
##     2  0 24  0  0  0
##     3  0  0 20  0  0
##     4  0  0  0 16  0
##     5  0  0  0  0  9
## [1] "Kappa overall =  1"
## [1] "Tau overall =  1"
## [1] "mean quality =  1"
## [1] "The quality is  1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect

## [1] "#########  Cramer's V =  1"
## [1] "classification error rate with altdata:  0"